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ANALYSIS  OF  LANDING-GEAR  BEHAVIOR  1 


By  Benjamin  Milwitzky  and  Francis  K.  Cook 


SUMMARY 

This  report  presents  a  theoretical  study  of  the  behavior  of  the 
conventional  type  of  oleo-pneumatic  landing  gear  during  the 
process  of  landing  impact.  The  basic  analysis  is  presented  in 
a  general  form  and  treats  the  motions  of  the  landing  gear  prior 
to  and  subsequent  to  the  beginning  of  shock-strut  deflection.  In 
the  analysis  of  the  first  phase  of  the  impact  the  landing  gear  is 
treated  as  a  si ngle-degree-of -freedom  system  in  order  to  deter¬ 
mine  the  conditions  of  motion  at  the  instan  t  of  in  itial  shock-strut 
deflection ,  after  which  instant  the  landing  gear  is  considered  as 
a  system  with  two  degrees  of  freedom.  The  equations  for  the 
two-degree-of -freedom  system  consider  such  factors  as  the 
hydraulic  (velocity  square )  resistance  of  the  orifice,  the  forces 
due  to  air  compression  and  internal  friction  in  the  shock  strut, 
the  nonlinear  force-deflection  characteristics  of  the  tire,  the  wing 
lift,  the  inclination  of  the  landing  gear,  and  the  effects  of  wheel 
spin-up  drag  loads. 

The  applicability  of  the  analysis  to  actual  landing  gears  has 
been  investigated  for  the  particular  case  of  a  vertical  landing 
gear  in  the  absence  of  drag  loads  by  comparing  calculated 
results  with  experimental  drop-test  data  for  impacts  with  and 
without  tire  bottoming.  The  calculated  behavior  of  the  landing 
gear  was  found  to  be  in  good  agreement  with  the  drop-test  data. 

Studies  have  also  been  made  to  determine  the  effects  of  varia¬ 
tions  in  such  parameters  as  the  dynamic  force-deflection 
characteristics  of  the  tire,  the  orifice  discharge  coefficient,  and  the. 
polytropic  exponent  for  the  air-compression  process,  which 
might  not  be  known  accurately  in  practical  design  problems. 

The  study  of  the  effects  of  variations  in  the  tire  characteristics 
indicates  that  in  the  case  of  a  normal  impact  without  tire, 
bottoming  reasonable  variations  in  the  force-deflection  character¬ 
istics  have  only  a  relatively  small  effect  on  the  calculated  behavior 
if  the  landing  gear.  Approximating  the  rather  complicated 
force-deflection  characteristics  if  the  actual  tire  by  simplified 
exponential  or  linear-segment  variations  appears  to  be  adequate 
for  practical  purposes.  Tire  hysteresis  was  found  to  be 
relatively  unimportant.  In  the  case  of  a  severe  impact  involving 
tire  bottoming,  flee  use  of  simplified  exponential  and  linear- 


segment  approximations  to  the  actual  tire  force-deflection 
characteristics,  which  neglect  the  effects  of  tire  bottoming, 
although  adequate  up  to  the  instant  of  bottoming,  fails  to  indicate 
the  pronounced  increase  in  landing-gear  load  that  results  from 
bottoming  of  the  tire.  The  use  of  exponential  and  linear- 
segment  approximations  to  the  tire  characteristics  which  take 
into  account  the  increased  stiffness  of  the  tire  which  results  from 
bottoming,  however,  yields  good  results. 

The  study  of  the  importance  of  the  discharge  coefficient  of  the 
orifice  indicates  that  the  magnitude  of  the  discharge  coefficient 
has  a  marked  effect  on  the  calculated  behavior  of  the  landing 
gear;  a  decrease  in  the  discharge  coefficient  (or  the  product  of  the 
discharge  coefficient  and  the  net  orifice  area)  results  in  an 
approximately  proportional  increase  in  the  maximum  upper- 
mass  acceleration. 

The  study  of  the  importance  of  the  air-compression  process 
in  the  shock  strut  indicates  that  the  air  springing  is  of  only 
minor  significance  throughout  most  of  the  impact  and  that 
variations  in  the  effective  polytropic  exponent  n  between  the 
isothermal  value,  of  1.0  and  the  near-adiabatic  value  of  1.3  hare 
only  a  secondary  effect  on  the  calculated  behavior  of  the  landing 
gear.  Even  the  assumption  of  constant  air  pressure  in  the  strut 
equal  to  the  initial  pressure,  that  is,  n—O,  yields  fairly  good 
results  which  may  he  adequate  for  many  practical  purposes. 

In  addition  to  the  more  exact  treatment,  an  investigation  has 
been  made  to  determine  the  extent  to  which  the -  basic  equations 
of  motion  can  be  simplified  and  still  yield  acceptable  results. 
This  study  indicates  that,  for  many  practical  purposes,  the 
air-pressure  force  in  the  shock  strut  can  V  completely  neglected, 
the  tire  force-deflection  relationship  ca  .  be  assumed  to  be  linear, 
and  the  lower  or  unsprung  mass  tan  be  taken  equal  to  zero. 
Generalization  of  the  equations  of  motion  for  this  simplified 
system  shows  that  the  behavior  of  the  system  is  completely 
determined  by  the  magnitude  of  one  parameter,  namely  the 
dimensionless  initial-velocity  parameter.  Solutions  of  these 
generalized  equations  ore  presented  in  terms  of  dimensionless 
variables  for  a  wide  range  of  landing-gear  and  impact  parameters 
which  may  be  useful  for  rapidly  estimating  lamling-gear 
performance  in  preliminary  design. 


1  Supersedes  NACA  TN  2755,  "Analysis  of  Landing-Gear  Behavior"  by  Benjamin  Milwitzky  and  Francis  E.  Took,  1032. 
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INTRODUCTION 

The  shock-absorbing  characteristics  of  airplane  landing 
gears  are  normally  developed  largely  by  means  of  extensive 
trial-and-error  drop  testing.  The  desire  to  reduce  the  ex¬ 
pense  and  time  required  by  such  methods,  as  well  as  to  pro¬ 
vide  a  more  rational  basis  for  the  prediction  of  wheel-inertia 
drag  loads  and  dynamic  stresses  in  flexible  airframes  during 
landing,  emphasizes  the  need  for  suitable  theoretical  methods 
for  the  analysis  of  landing-gear  behavior.  Such  theoretical 
methods  should  And  application  in  the  design  of  landing 
gears  and  complete  airplane  structures  by  permitting 

(a)  the  determination  of  the  behavior  of  a  given  landing- 
gear  configuration  under  varying  impact  conditions  (velocity 
at  contact,  weight,  wing  lift,  etc.) 

(b)  the  development  of  a  landing-gear  configuration  to 
obtain  a  specified  behavior  under  given  impact  conditions 

(c)  a  more  rational  approach  to  the  determination  of  wheel 
spin-up  and  spring-back  loads  which  takes  into  account  the 
shock-absorbing  characteristics  of  the  particular  landing  gear 
under  consideration 

(d)  improved  determination  of  dynamic  loads  in  flexible 
airplane  structures  during  landing.  This  problem  may  be 
treated  either  by  calculating  the  response  of  the  elastic  sys¬ 
tem  to  landing-gear  forcing  functions  determined  under  the 
assumption  that  the  airplane  is  a  rigid  body  or  by  the  simul¬ 
taneous  solution  of  the  equations  of  motion  for  the  landing 
gear  coupled  with  the  equations  representing  the  additional 
degrees  of  freedom  of  the  structure.  In  many  cases  the  former 
approach  should  be  sufficiently  accurate,  but  in  some 
instances,  particularly  when  the  landing-gear  attachment 
points  experience  large  displacements  relative  to  the  nodal 
points  of  the  flexible  system,  the  latter  approach,  which  takes 
into  account  the  interaction  between  the  deformations  of  the 
structure  and  the  landing  gear,  may  be  required  in  order  to 
represent  the  system  adequately. 

Since  many  aspects  of  the  landing-impact  problem  are  so 
intimately  connected  with  the  mechanics  of  the  landing  gear, 
the  subject  of  landing-gear  behavior  has  received  analytical 
treatment  at  various  times  (see  bibliography).  Many  of 
the  earlier  investigations,  in  order  to  reduce  the  mathematical 
complexity  of  the  analysis,  were  limited  to  consideration  of 
highly  simplified  linear  systems  which  have  little  relation  to 
practical  landing  gears.  Some  of  the  more  recent  papers 
consider,  with  different  degrees  of  simplification,  more  real¬ 
istic  nonlinear  systems.  The  present  report  represents  an 
attempt  at  a  more  complete  analysis  of  the  mechanics  of 
practical  landing  gears  and,  in  addition,  investigates  the  im¬ 
portance  of  the  various  elements  which  make  up  the  landing 
gear,  as  well  as  the  extent  to  which  the  system  can  be  reason¬ 
ably  simplified  for  the  purpose  of  rapid  analysis. 

The  basic  analysis  is  presented  in  a  general  form  and  takes 
into  account  such  factors  as  the  hydraulic  (velocity  square) 
resistance  of  the  orifice,  the  forces  due  to  air  compression  and 
internal  friction  in  the  shock  strut,  the  nonlinear  force- 
deflection  characteristics  of  the  tire,  the  wing  lift,  the  inclina¬ 
tion  of  the  landing  gear  and  the  effects  of  wheel  spin-up  drag 
loads.  An  evaluation  of  the  applicability  of  the  analysis  to 
actual  landing  gears  is  presented  for  the  case  of  a  vertical 
landing  gear  in  the  absence  of  drag  loads  by  comparing  cal¬ 
culated  results  with  drop-test  data. 


Since  some  parameters,  such  as  the  dynamic  force- 
deflection  characteristics  of  the  tire,  the  orifice  discharge 
coefficient,  and  the  poly  tropic  exponent  for  the  air-compression 
process,  may  not  be  accurately  known  in  practical  design 
problems,  a  study  is  made  to  assess  the  effects  of  variations 
in  these  parameters  on  the  calculated  landing-gear  behavior. 

Studies  are  also  presented  to  evaluate  the  extent  to  which 
the  dynamical  system  can  be  simplified  without  greatly  im¬ 
pairing  the  validity  of  the  calculated  results.  In  addition  to 
the  investigations  for  specific  cases,  generalized  solutions  for 
the  behavior  of  a  simplified  system  are  presented  for  a  wide 
range  of  landing-gear  and  impact  parameters  which  may  be 
useful  in  preliminary  design. 

SYMBOLS 

Aa  pneumatic  area 

At,  hydraulic  area 

A„  area  of  opening  in  orifice  plate 

Ai  internal  cross-sectional  area  of  shock-strut  inner 

cylinder 

At  external  cross-sectional  area  of  shock-strut  inner 

cylinder 

Ap  cross-sectional  area  of  metering  pin  or  rod  in 

plane  of  orifice 
An  net  orifice  area 

Cd  orifice  discharge  coefficient 

d  overall  diameter  of  tire 

Fa  pneumatic  force  in  shock  strut 

F„  hydraulic  force  in  shock  strut 

Fr  friction  force  in  shock  strut 

Fs  total  axial  shock-strut  force 

Ft  normal  force  on  upper  bearing  (attached  to  inner 

cylinder) 

F2  normal  force  on  lower  bearing  (attached  to  outer 

cylinder) 

Fya  force  normal  to  axis  of  shock  strut,  applied  at 

axle 

FVa  vertical  force,  applied  at  axle 

F„a  horizontal  force,  applied  at  axle 

FRa  resultant  force,  applied  at  axle 

FSf  force  parallel  to  axis  of  shock  strut,  applied  to 

tire  at  ground 

F»  force  normal  to  axis  of  shock  strut,  applied  to 

tire  at  ground 

FVf  vertical  force,  applied  to  tire  at  ground 

F„t  horizontal  force,  applied  to  tire  at  ground 

Fr  resultant  force,  applied  to  tire  at  ground 

g  gravitational  constant 

Kl  lift  factor,  L/W 

L  lift  force 

l,  axial  distance  between  upper  and  lower  bearings, 

for  fully  extended  shock  strut 

It  axial  distance  between  axle  and  lower  bearing 

(attached  to  outer  cylinder),  for  fully  ex¬ 
tended  shock  strut 

a,b,m,r  constants  corresponding  to  the  various  regimes 
of  the  tire-deflection  process 
a'  combined  constant,  ad 

m'  combined  constant,  mdr 

n  poly  tropic  exponent  for  air- compression  process 

in  shock  stint 
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li  Reynolds  number 

pa  air  pressure  in  upper  chamber  of  shock  strut 

pk  hydraulic  pressure  in  lower  chamber  of  shock 

strut 

Q  volumetric  rate  of  discharge  through  orifice 

rd  radius  of  deflected  tire 

s  shock-strut  axial  stroke 

T  wheel  inertia  torque  reaction 

t  time  after  contact 

r  time  after  beginning  of  shock-strut  deflection 

v  air  volume  of  shock  strut 

Iw  polar  moment  of  inertia  for  wheel  assembly 

about  axle 

T’v-  vertical  velocity 

TV  horizontal  velocity 

11 '  total  dropping  weight 

H’,  weight  of  upper  mass  above  strut 

li’.  weight  of  lower  mass  below  strut 

it  horizontal  displacement  of  lower  mass  from 

position  at  initial  contact 
vertical  displacement  of  upper  mass  from  posi¬ 
tion  at  initial  contact 

z«  vertical  displacement  of  lower  mass  from  posi¬ 

tion  at  initial  contact 

iii  dimensionless  upper-mass  displacement  from 

position  at  initial  contact 

iii  dimensionless  lower-mass  displacement  from 

position  at  initial  contact 

<r  dimensionless  shock-strut  stroke,  V\-ut 

0  dimensionless  time  after  contact 

y  angle  between  shock-strut  axis  and  vertical 

I  Ui"  do 

rj,  shock-strut  effectiveness,  - - 

“1  maximal 


Vie 


M 

Ml 


M2 

P 

a 

Axes: 

z 

X 

Subscripts: 

0 

T 

su 

mar 

Notation: 


landing-gear  effect  iveness, 


/: 


Ml 


"  dUi 


time  interval  in  numerical  integration  procedures 
coefficient  of  friction  between  tire  and  runway 
coefficient  of  friction  for  upper  bearing  (attached 
to  inner  cylinder) 

coefficient  of  ftiction  for  lower  bearing  (attached 
to  outer  cylinder) 
mass  density  of  hydraulic  fluid 
angular  acceleration  of  wheel 


vertical  axis,  positive  downward 
horizontal  axis,  positive  rearward 


at  instant  of  initial  contac  t 
at  instant  of  initial  shock-strut  deflection 
at  instant  of  wheel  spin-up 
maximum  value 


|(  )|  absolute  value  of  (  ) 

(  )*  estimated  value  of  (  ) 

The  use  of  dots  over  symbols  indicates  differentiation  with 
respect  to  time  l  or  t! 

Prime  marks  indicate  differentiation  with  respect  to 
dimensionless  time  8. 


MECHANICS  OF  LANDING  GEAR 

DYNAMICS  OF  SYSTEM 

In  view  of  the  fact  that  landing-gear  performance  appears 
to  be  relatively  unaffected  bv  the  elastic  deformations  of 
the  airplane  structure  (see,  for  example,  refs.  1  and  2)  par¬ 
ticularly  since  in  many  cases  the  main  gears  arc*  located 
fairly  close  to  the  nodal  points  of  the  fundamental  bending 
mode  of  the  wing,  that  part  of  the  airplane  which  acts  on  a 
given  gear  can  generally  be  considered  as  a  rigid  mass. 
As  a  result,  landing-gear  drop  tests  are  often  conducted  in 
a  jig  where  the  mass  of  the  airplane  is  represented  by  a 
concentrated  weight.  In  particular  instances,  however,  such 
as  in  the  case  of  airplanes  having  large  concentrated  masses 
disposed  in  an  outboard  position  in  the  wings,  especially 
airplanes  equipped  with  bicycle  landing  gear,  consideration 
of  the  interaction  between  the  deformation  of  the  airplane 
structure  and  the  landing  gear  may  be  necessary  to  repre¬ 
sent  the  system  adequately. 

Since  the  present  report  is  concerned  primarily  with  the 
mechanics  of  the  landing  gear,  it  is  assumed  in  the  analysis 
that  the  landing  gear  is  attached  to  a  rigid  mass  which  has 
freedom  only  in  vertical  translation.  The  gear  is  assumed 
infinitely  rigid  in  bending.  The  combination  of  airplane 
and  landing  gear  considered  therefore  constitutes  a  system 
having  two  degrees  of  freedom  (see  fig.  1(a))  as  defined  by 
the  vertical  displacement  of  the  upper  mass  and  the  vertical 
displacement  of  the  lower  or  unsprung  mass,  which  is  also 
the  tire  deflection.  The  strut  stroke  .«  is  determined  by 
the  difference  between  the  displacements  z ,  and  z2  and.  in 
the  case  of  inclined  gears,  by  the  angle  y  between  the  axis 
of  the  strut  and  the  vertical.  For  inclined  gears,  compression 
of  the  shock  strut  produces  a  horizontal  displacement  of  the 
axle  j'2.  From  consideration  of  the  kinematics  of  the  system 

it  can  be  seen  that  s—~-  —  and  jy=.v  sin  y  =  (V  —  s2)tan  y. 
cos  y 

In  the  analysis,  external  lift  forces,  corresponding  to  the 
aerodynamic  lift,  are  assumed  to  act  on  the  system  through¬ 
out  the  impact.  In  addition  to  the  vertical  forces,  arbitrary 
drag  loads  are  considered  to  act  between  the  tire  and  the 
ground. 

The  system  treated  in  the  analysis  may  therefore  be  con¬ 
sidered  to  represent  either  a  landing-gear  drop  test  in  a  jig 
where  wing  lift  and  drag  loads  are  simulated,  or  the  landing 
impact  of  a  rigid  airplane  if  rotational  motions  are  neglected. 
Rotational  freedom  of  the  airplane,  where  significant,  may 
be  taken  into  account  approximately  by  use  of  an  appro¬ 
priate  effective  mass  in  the  analysis. 

Figure  1  (b)  shows  a  schematic  representation  of  a  typical 
oleo-pneumatic  shock  strut  used  in  American  practice.  The 
lower  chamber  of  the  strut  contains  hydraulic  fluid  and  the 
upper  chamber  contains  air  under  pressure.  The  outer  cyl¬ 
inder  of  the  strut,  which  is  attached  to  the  upper  mass, 
contains  a  perforated  tube  which  supports  a  plate  with  a 
small  orifice,  through  which  the  hydraulic,  fluid  is  forced  to 
flow  at  high  velocity  as  a  result  of  the  telescoping  of  the 
strut.  The  hydraulic  pressure  drop  across  the  orifice  thus 
produced  resists  the  closure  of  the  strut,  ami  the  turbulence 
created  provides  a  powerful  means  of  absorbing  and  dis¬ 
sipating  a  large  part  of  the  impact  energy.  In  some  struts 
the  orifice  area  is  constant ;  whereas,  in  other  cases  a  metering 
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(a)  System  with  two  degrees  of  freedom.  (,!>)  Schematic  representation  of  shock  strut. 

Fioi  rk  1. —  Dynamical  system  considered  in  analysis. 


pin  or  rod  is  used  to  control  the  size  of  the  orifice  and  govern 
the  performance  of  the  strut. 

The  compression  of  the  strut  produces  an  increase  in  the 
air  pressure  which  also  resists  the  closure  of  the  strut.  In 
figure  1  (b)  p„  represents  the  oil  pressure  in  the  lower  chamber 
and  pa  represents  the  air  pressure  in  the  upper  chamber. 

In  addition  to  the  hydraulic  resistance  and  air-pressure 
foices,  internal  bearing  friction  also  contributes  forces  which 
can  appreciably  affect  the  behavior  of  the  strut. 

The  forces  created  within  the  strut  impart  an  acceleration 
to  the  upper  mass  and  also  produce  an  aeeeleration  of  the 
lower  mass  and  a  deflection  of  the  tire.  Figure  1  (c)  shows 
the  balance  of  forces  ami  reactions  for  the  wheel,  the  inner 
cylinder,  and  the  outer  cylinder.  It  is  clear  that  the  strut 


and  the  tire  mutually  influence  tin1  behavior  of  one  another 
and  must  be  considered  simultaneously  in  analyzing  t  he  syst  <•  m. 

FORCES  IN  SHOCK  STRUT 

From  consideration  of  the  pressures  acting  in  the  shock 
strut  it  can  be  readily  seen  from  figure  1  (b)  that  the  total 
axial  force  due  to  hydraulic  resistance,  air  compression,  and 
hearing  friction  can  be  expressed  by 

F.s  P, ,(zli  -1  ,J )  "I  Paf-'ls  -'ll)  “h  Pa- 1,, T  * 

where 

-4 1  internal  cross-sectional  area  of  inner  cylinder 

-12  external  cross-sectional  area  of  inner  cylinder 

A„  cross-sectional  area  of  metering  pin  or  rod  in  plane 
of  orifice 
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This  expression  can  also  be  written  as 

f's—(Ph — Pa){A\— Ar  paA«-\-  F j 
=  (Ph—  Pa)A.h+J>aAa+Ff 
=  Fk+F.+F,  (1) 

where 

Ph—Pa  pressure  drop  across  the  orifice 

A/,  hydraulic  area  (A,—Ap  for  the  strut  shown  in  fig.  1) 

A„  pneumatic  area  {At  for  the  strut  shown  in  fig.  1 ) 

In  this  report  the  terms  (p*—p*)A*  and  pdAd  are  referred 
to  as  hydraulic  force  Fh  and  pneumatic  force  Fa,  respec¬ 
tively.  For  the  strut  shown  in  figure  1,  the  hydraulic 
and  pneumatic  areas  are  related  to  the  strut  dimensions  as 
previously  noted.  In  the  case  of  struts  having  different 
internal  configurations,  the  hydraulic  and  pneumatic  areas 
may  bear  somewhat  different  relations  to  the  dimensions 
of  the  strut.  In  such  cases,  however,  consideration  of  the 
pressures  acting  on  the  various  components  of  the  strut 
should  permit  these  areas  to  be  readily  defined. 

Hydraulic  force. — The  hydraulic  resistance  in  the  shock 
strut  results  from  the  pressure  difference  associated  with  the 
flow  through  the  orifice.  In  a  landing  gear  the  orifice  area 
is  usually  small  enough  in  relation  to  the  diameter  of  the 
strut  so  that  the  jet  velocities  and  Reynolds  numbers  are 
sufficiently  large  that  the  flow  is  fully  turbulent.  As  a 
result  the  damping  force  varies  as  the  square  of  the  tele¬ 
scoping  velocity  rather  than  linearly  with  the  velocity. 
Since  the  hydraulic  resistance  is  the  major  component  of 
the  total  shock-strut  force,  viscous  damping  cannot  be 
reasonably  assumed,  even  though  such  an  assumption 
would  greatly  simplify  the  analysis. 

The  hydraulic  resistance  can  be  readily  derived  by  making 
use  of  the  well-known  equation  for  the  discharge  through 
an  orifice,  namely, 

Q  =  C„An^\  (p»-~pa) 

where 

0  volumetric  rate  of  discharge 
C„  coefficient  of  discharge 
An  net  orifice  area 

ph  hydraulic  pressure  in  lower  chamber 
p„  air  pressure  in  upper  chamber 
P  mass  density  of  hydraulic  fluid 
From  considerations  of  continuity,  the  volumetric  rate  of 
discharge  can  also  be  expressed  as  the  product  of  the  tele¬ 
scoping  velocity  8  and  the  hydraulic  area  Ah 

Q  —  A/,'* 

Equating  the  preceding  expressions  for  the  discharge  per¬ 
mits  writing  the  following  simple  equation  for  the  pressure 
drop  across  the  orifice 

_  pAh2s2 

Ph  Pa~2  {C„Any 

The  hydraulic  resistance  Fh  due  to  the  telescoping  of  the 
strut  is  given  by  the  product  of  the  differential  pressure 


ph—pa  and  the  area  .1*  which  is  subjected  to  the  hydraulic 
pressure,  as  previously  noted.  Thus 


F_  PA„3 
*  2{CdAny 


(2) 


Equation  (2)  can  be  made  applicable  to  both  the  compres¬ 
sion  and  elongation  strokes  by  introducing  the  factor  rr-v 
to  indicate  the  sign  of  the  hydraulic  resistance ;  thus 


l~i _  £  P-Aa3^ 

2{CdAnY  h 


(2a) 


The  net  orifice  area  A*  may  be  either  a  constant  or,  when  a 
metering  pin  is  used,  can  vary  with  strut  stroke;  that  is, 
An—A0—Ap=A„(s),  where  A0  is  the  area  of  the  opening  in 
the  orifice  plate  and  Ap  is  the  area  of  the  metering  pin  in  the 
plane  of  the  orifice.  At  the  present  time  there  appears  to  be 
some  tendency  to  eliminate  the  metering  pin  and  use  a  con¬ 
stant  orifice  area,  particularly  for  large  airplanes,  in  which 
case  An=A„.  In  the  general  case,  the  orifice  discharge 
coefficient  might  be  expected  to  vary  somewhat  during  an 
impact  because  of  changes  in  the  size  and  configuration  of 
the  net  orifice  area,  changes  in  the  exit  conditions  on  the 
downstream  face  of  the  orifice  due  to  variations  in  the  amount 
of  hydraulic  fluid  above  the  orifice  plate,  changes  in  the  entry 
conditions  due  to  variations  in  the  length  of  the  flow  chamber 
upstream  of  the  orifice,  and  because  of  variations  in  the 
Reynolds  number  of  the  flow,  so  that,  in  general,  Cd  —  Cd(n,  R). 
Although  the  individual  effects  of  these  factors  on  the  dis¬ 
charge  coefficients  for  orifices  in  shock  struts  have  not  been 
evaluated,  there  is  some  experimental  evidence  to  indicate 
appreciable  variations  of  the  discharge1  coefficient  during 
impact,  particularly  in  the  case  of  struts  with  metering  pins. 
It  might  be  expected  that  such  variations  would  be  con¬ 
siderably  smaller  for  gears  having  a  constant  orifice  area. 

In  order  to  evaluate  the  precision  with  which  the  orifice 
discharge  coefficient  has  to  be  known,  a  brief  study  is 
presented  in  a  subsequent  section  which  shows  the  effect  of 
the  discharge  coefficient  on  the  calculated  behavior  of  a 
landing  gear  with  a  constant  orifice  area,  under  the  assump¬ 
tion  that  the  discharge  coefficient  is  constant  during  the 
impact . 

The  foregoing  discussion  has  been  concerned  primarily  with 
the  compression  stroke  of  the  shock  strut.  Most  struts 
incorporate  some  form  of  pressure-operated  rebound  check 
valve,  sometimes  called  a  snubber  valve,  which  comes  into 
action  after  the  maximum  stroke  has  been  attained  and  closes 
off  the  main  orifice  as  soon  as  the  strut  begins  to  elongate,  so 
that  the  fluid  is  forced  to  return  to  the  lower  chamber  through 
small  passages.  The  action  of  the  snubber  valve  introduces 
greatly  increased  hydraulic  resistance  to  dissipate  the  energy 
stored  in  the  strut  in  the  form  of  air  pressure  and  to  prevent 
excessive  rebound.  The  product  CdAn  to  be  used  in  equation 
(2a)  during  the  elongation  stroke  is  generally  uncertain.  The 
exact  area  An  during  elongation  is  usually  somewhat  difficult 
to  define  from  the  geometry  of  the  strut  since  in  many  cases 
the  number  of  connecting  passages  varies  w  ith  stroke  and  the 
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leakage  area  around  the  piston  may  be  of  the  same  order  of 
magnitude  as  the  area  of  the  return  passages.  Furthermore, 
the  magnitude  of  the  orifice  discharge  coefficient,  and  even 
possibly  the  nature  of  the  resistance,  are  questionable  due  to 
the  foaming  state  of  the  returning  fluid.  Fortunately,  the 
primary  interest  is  in  the  compression  process  rather  than 
the  elongation  process  since  the  maximum  load  always  occurs 
before  the  maximum  strut  stroke  is  reached. 

Pneumatic  force. — The  air-pressure  force  in  the  upper 
chamber  is  determined  by  the  initial  strut  inflation  pres¬ 
sure,  the  area  subjected  to  the  air  pressure  (pneumatic  area), 
and  the  instantaneous  compression  ratio  in  accordance 
with  the  polytropic  law  for  compression  of  gases,  namely 
p„vn= Constant,  or 


Pa=P% 


n 


where 

pa  air  pressure  in  upper  chamber  of  shock  strut 

pa o  air  pressure  in  upper  chamber  for  fully  extended  strut 

v  air  volume  of  shock  strut 

r0  air  volume  for  fully  extended  strut 

Since  the  instantaneous  air  volume  is  equal  to  the  difference 

between  the  initial  air  volume  and  the  product  of  the  stroke 


and  pneumatic  area  Aa, 


P“_P“°G0  -A.*) 


The  force  due 


to  the  air  pressure  is  simply  the  product  of  the  pressure  and 
the  pneumatic  area: 


Fa-p^Aa  (r#_I^)  (3) 


In  the  preceding  equations,  the  effective  polytropic 
exponent  n  depends  on  the  rate,  of  compression  and  the  rate 
of  heat  transfer  from  the  air  to  the  surrounding  environment. 
Low  rates  of  compression  would  be  expected  to  result  in 
values  of  n  approaching  the  isothermal  value  of  1.0;  whereas 
higher  values  of  n,  limited  by  the  adiabatic  value  of  1.4, 
would  be  expected  for  higher  rates  of  compression.  The 
actual  thermodynamic  process  is  complicated  by  the  violent 
mixing  of  the  highly  turbulent  efflux  of  hydraulic  fluid  and 
the  air  in  the  upper  chamber  during  impact.  On  the  one 
hand,  the  dissipation  of  energy  in  the  production  of  turbu¬ 
lence  generates  heat;  on  the  other  hand,  heat  is  absorbed  by 
the  aeration  and  vaporization  of  the  fluid.  The  effect  of  this 
mixing  phenomenon  on  the  polytropic  exponent  or  on  the 
equivalent  air  volume  is  not  clear.  A  limited  amount  of 
experimental  data  obtained  in  drop  tests  (refs,  :i  and  4), 
however,  indicates  that  the  effective  polytropic  exponent 
may  be  in  the  neighborhood  of  1.1  for  practical  cases.  A 
brief  study  of  the  importance  of  the  air-compression  process 
and  the  effects  which  different  values  of  n  may  have  on  the 
calculated  behavior  of  the  landing  gear  is  presented  in  a 
subsequent  section. 

Internal  friction  force. — In  the  literature  on  machine  de¬ 
sign  the  wide  range  of  conditions  under  which  frictional 
resistance  can  occur  between  sliding  surfaces  is  generally 
classified  in  three  major  categories,  namely,  friction  betwe  -n 
dry  surfaces,  friction  between  imperfectly  lubricated  surfaces, 
and  friction  between  perfectly  lubricated  surfaces.  In  the 


case  of  dry  friction,  the  resistance  depends  on  the  physical 
characteristics  of  the  sliding  surfaces,  is  essentially  propor¬ 
tional  to  the  normal  force,  and  is  approximately  independent 
of  the  surface  area.  The  coefficient  of  friction  n,  defined  as 
the  ratio  of  the  frictional  resistance  to  the  normal  force,  is 
generally  somewhat  greater  under  conditions  of  rest  (static 
friction)  than  under  conditions  of  sliding  (kinetic  friction). 
Although  the  coefficient  of  kinetic  friction  generally  de¬ 
creases  slightly  with  increasing  velocity,  it  is  usually  con¬ 
sidered,  in  first  approximation,  to  be  independent  of  velocity. 
If,  on  the  other  hand,  the  surfaces  are  completely  separated 
by  a  fluid  film  of  lubricant,  perfect  lubrication  is  said  to  exist. 
Under  these  conditions  the  resistance  to  relative  motion 
depends  primarily  on  the  magnitude  of  the  relative  velocity, 
the  physical  characteristics  of  the  lubricant,  the  area,  and 
the  film  thickness,  and  is  essentially  independent  of  the 
normal  force  and  the  characteristics  of  the  sliding  surfaces. 
Perfect  lubrication  is  rarely  found  in  practice  but  is  most 
likely  under  conditions  of  high  velocity  and  relatively  small 
normal  pressure,  where  the  shape  of  the  sliding  surfaces  is 
conducive  to  the  generation  of  fluid  pressure  by  hydro- 
dynamic  action.  In  most  practical  applications  involving 
lubrication,  a  state  of  imperfect  lubrication  exists  and  the 
resistance  phenomenon  is  intermediate  between  that  of  dry 
friction  and  perfect  lubrication. 

In  the  case  of  landing-gear  shock  struts,  the  conditions 
under  which  internal  friction  is  of  concern  usually  involve 
relatively  high  normal  pressures  and  relatively  small  sliding 
velocities.  Moreover,  the  usual  types  of  hydraulic  fluid 
used  in  shock  struts  have  rather  poor  lubricating  properties, 
and  the  shape  of  the  bearing  surfaces  is  generally  not  con¬ 
ducive  to  the  generation'of  hydrodynamic  pressures.  It 
would  therefore  appear  that  the  lubrication  of  shock  strut 
bearings  is,  at  best,  imperfect;  in  fact,  the  conditions  appear 
to  approach  closely  those  for  dry  friction.  In  the  present 
analysis,  therefore,  it  is  assumed,  in  first  approximation,  that 
the  internal  friction  between  the  bearings  and  the  cylinder 
walls  follows  laws  similar  to  those  for  dry  friction;  that  is, 
the  friction  force  is  given  by  the  product  of  the  normal  force 
and  a  suitably  chosen  coefficient  of  friction. 

With  these  assumptions  the  internal  friction  forces  pro¬ 
duced  in  the  strut  depend  on  the  magnitude  of  the  forces  on 
the  axle,  the  inclination  of  the  gear,  the  spacing  of  the  bear¬ 
ings,  and  the  coefficient  of  friction  between  the  bearings  and 
the  cylinder  walls.  Figure  1(c)  schematically  illustrates  the 
balance  of  forces  acting  on  the  various  components  of  the 
landing  gear.  The  total  axial  friction  in  the  shock  strut  is 
the  sum  of  the  friction  forces  contributed  by  each  of  the 
bearings: 

F/-||y(M,|F1|  +  M.|i',.|) 

where 

F,  axial  friction  force 

lii  coefficient  of  friction  for  upper  bearing  (attached  to 
inner  cylinder) 

F,  normal  force  on  upper  bearing  (attached  to  inner 
cylinder) 

to  coefficient  of  friction  for  lower  bearing  (attached  to 
outer  cylinder) 
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normal  force  on  lower  bearing  (attached  to  outer 
cylinder) 

factor  to  indicate  sign  of  friction  force 


During  the  interval  prior  to  the  beginning  of  shock-strut 
motion  the  friction  forces  depend  on  the  coefficients  of  static 
friction;  after  the  strut  begins  to  telescope  the  coefficients  of 
kinetic  friction  apply. 

From  considerations  of  the  balance  of  moments  it  can  be 
seen  from  figure  1  (c)  that 


so  that 

F/^W\  l  /r+s+fl2]  (4) 

where 

F.va=  FVg  sin  F„a  cos  f  (4a) 

and 

FSa  force  normal  tc  strut  applied  at  axle 
F\a  vertical  force  applied  at  axle 
F„a  horizontal  force  applied  at  axle 
f  angle  between  strut  axis  and  vertical 
/,  axial  distance  between  upper  and  lower  bearings,  for 
fully  extended  strut 

/•>  axial  distance  between  axle  and  lower  bearing  (attached 
to  outer  cylinder),  for  fully  extended  strut 
The  quantities  FXa,  FVa,  and  F„a  are  forces  applied  at  the 
axle  and  differ  from  the  ground  reactions  by  amounts  equal 
to  the  inertia  forces  corresponding*to  the  respective  accelera¬ 
tion  component  of  the  lower  mass.  Since  the  inner  cylinder 
generally  represents  only  a  relatively  small  fraction  of  the 
lower  mass,  the  lower  mass  may  reasonably  be  assumed  to  lx 
concentrated  at  the  axle.  With  this  assumption,  the  rela¬ 
tionships  between  the  forces  at  the  axle  and  the  forces  at  the 
ground  are  given  by 


The  normal  force  at  the  axle  can  therefore  be  expressed  in 
terms  of  the  ground  reactions  and  the  component  accelera¬ 
tions  of  the  lower  mass  by 


Fya—(^  Fvt-\-  ' ‘  z-i—  ttV)  sin  f  —  (  F„  —  '  j;2)  cos  (4b) 

where 

Fv  vertical  force  applied  to  tire  at  ground 
F„g  horizontal  force  applied  to  tire  at  ground 

effective  mass  below  shock  strut,  assumed  concent  rated 
g  at  axle 

x«  horizontal  acceleration  of  axle 

Z'i  vertical  acceleration  of  axle 

In  the  case  of  an  inclined  landing  gear  having  infinite  stiff¬ 
ness  in  bending,  the  horizontal  displacement  of  the  lower  mass 
i2  is  related  to  the  vertical  displacements  of  the  upperand  lower 


masses  by  the  kinematic  relationship  x-,—  (a  —  5>)tan  <p,  as 
previously  noted.  Double  differentiation  of  this  relation¬ 
ship  gives  x2=  (ij, —  z2)  tan  if.  Substitution  of  this  expression 
into  equation  (4b)  gives 


FN  —  Fv.  sin  if—  Fv.  cos  <p- 1 - -  z,  sin  <f—  U’2  sin  <f  (4c) 

“  •  I  ^ 


In  equation  (4c)  (he  quantity  zx  sin  <f  represents  the  ac¬ 
celeration  of  the  lower  mass  normal  to  the  strut  axis  when  the 
gear  is  rigid  in  bending.  In  the  case  of  a  gear  flexible  in 
bending,  the  normal  acceleration  of  the  lower  mass  is  not 
eompletely  determined  bv  the  vertical  acceleration  of  the 
upper  mass  and  the  angle  of  inclination  of  the  gear.  If  it 
should  be  necessary  to  take  into  account,  in  particular  cases, 
the  effects  of  gear  flexibility  on  the  relationship  between  the 
normal  force  on  the  axle  and  the  ground  reactions,  the  quan¬ 
tity  z i  sin  f  in  equation  (4c)  may  be  replaced  by  estimated 
values  of  the  actual  normal  acceleration  of  the  lower  mass  as 
determined  from  consideration  of  the  bending  response  of  the 
gear  to  the  applied  force's  normal  to  the  gear  axis.  The  effects 
of  gear  flexibility  are  not  considered  in  more  detail  in  the 
present  analysis. 


FORCES  ON  TIRE 

Figure  2  (a)  shows  dynamic  force-deflection  characteris¬ 
tics  for  a  27-inch  smooth-contour  (type  I)  tire  inflated  to  32 
pounds  per  square  inch.  These  characteristics  were  deter¬ 
mined  from  time-history  measurements  of  vertical  ground 
force  and  tire  deflection  in  landing-gear  drop  tests  with  a 
nonrotating  wheel  at  several  vertical  velocities.  As  can  be 
seen,  the  tire  compresses  along  one  curve  and  unloads  along 
another,  the  hysteresis  loop  indicating  appreciable  energy 
dissipation  in  the  tire.  There  is  some  question  as  to  whether 
the  amount  of  hysteresis  would  be  as  great  if  the  tire  were 
rotating,  as  in  a  landing  with  forward  speed.  The  force- 
deflection  curve  for  a  velocity  of  1 1.63  feet  per  second  is  for 
a  severe  impact  in  which  tire  bottoming  occurs  and  shows 
the  sharp  increase  in  force  with  deflection  subsequent  to 
bottoming. 

In  figure  2  (b)  the  same  force-deflection  characteristics 
are  shown  plotted  on  logarithmic  coordinates.  As  can  be 
seen,  the  force  exhibits  an  exponential  variation  with  deflec¬ 
tion.  A  systematized  representation  of  the  force-deflection 
relationship  can  therefore  be  obtained  by  means  of  simple 
equations  having  the  form 


where 

Fy  vertical  force,  applied  to  tire  at  ground 

vertical  displacement  of  lower  mass  from  position  at 
initial  contact  (radial  deflection  of  tire) 

(l  overall  diameter  of  tire 

m,  r  constants  corresponding  to  the  various  regimes  of  the 
tire-deflection  process 
m'  combined  constant,  mil' 
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(a)  I'niform  coordinate*.  (b)  Logarithmic  coordinates’. 

Fiovre  2. — Dynamic  force-deflection  cliaracteristics  of  tire. 


It  may  be  noted  from  figure  2  that  essentially  the  same 
force-deflection  curve  holds  during  compression  for  all  impact 
velocities,  up  to  the  occurrence  of  tire  bottoming,  and  that  in 
figure  2  (b)  the  slopes  of  the  curves  in  each  of  the  several 
regimes  of  the  tire-deflection  process  are  also  independent  of 
velocity,  except  in  the  compression  regime  following  tire 
bottoming. 

Figure  2  also  shows  simple  approximations  to  the  tire 
characteristics  which  were  obtained  by  fitting  straight-line 
segments  (long-dashed  lines)  to  the  actual  force-deflection 
curves  in  figure  2  (a)  for  impacts  at  8. SO  and  11.63  feet  pet- 
second.  These  approximations,  hereinafter  referred  to  as 
linear-segment  approximations,  are  included  in  a  study, 
presented  in  a  subsequent  section,  to  evaluate  the  degree  of 
accuracy  required  for  adequate  representation  of  the  tire 
characteristics.  The  various  representations  of  the  tire 
characteristics  considered  and  the  pertinent  constants  for 
each  regime  of  tire  deflection  are  shown  in  figure  3. 

EQUATIONS  OF  MOTION 

'Fhe  internal  axial  force  Fs  produced  by  the  shock  strut 
was  shown  in  a  previous  section  to  be  equal  to  the  sum  of  the 


hydraulic,  pneumatic,  and  friction  forces,  as  given  In¬ 
equation  (1).  Since  these  forces  act  along  the  axis  of  the 
strut,  which  may  be  inclined  to  the  vertical  by  an  angle  <p, 
the  vertical  component  of  the  axial  shock-strut  force  is  given 
In  Fs  cos  <p.  The  vertical  component  of  the  force  normal 
to  the  shock  strut  is  given  by  FK  sin  <p.  These  forces  act 
in  conjunction  with  the  lift  force  and  weight  to  produce  an 
acceleration  of  the  upper  mass.  The  equation  of  motion 
for  the  upper  mass  is 


Fs  cos  g>+  /-.va  sin  *p-\-L  —  II  i  — 


(6) 


The  vertical  components  of  the  axial  and  normal  shock- 
strut  forces  also  act,  in  conjunction  with  the  weight  of  the 
lower  mass,  to  produce  a  deformation  of  the  tire  and  an 
acceleration  of  the  lower  mass.  The  equation  of  motion 
for  the  lower  mass  is 


Fs  cos  <p  +  Fs.  sin  <p  +  M  2 


W2 


where  the  vertical  ground  reaction  F 


z2 — Fvt(z2)  (7) 

r  is  expressed  as  a 
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Exponential  variation; 


I3xi0j 


/ 


(a)  Impact  without  tire  bottoming,  IV0=8.8(i  feet  per  second.  (b)  Impact  with  tire  bottoming,  IV#=  11.63  feet  per  second. 

Figure  3. — Tire  characteristics  considered  in  solutions  (logarithmic  coordinates). 


function  of  the  tire  deflection  z2.  The  relationship  between 
Fvt  and  z2  has  been  discussed  in  the  previous  section  on  tire 
characteristics. 

By  combining  equations  (6)  and  (7).  the  vertical  ground 
force  can  be  written  in  terms  of  the  inertia  reactions  of  the 
upper  and  lower  masses,  the  lift  force,  and  the  total  weight. 
The  overall  dynamic  equilibrium  is  given  by 

Fv,W=-y5,-yM+F  (8) 

MOTION  PRIOR  TO  SHOCK-STRUT  DEFLECTION 

Conventional  oleo-pneumatic  shock  struts  are  inflated 
to  some  finite  pressure  in  the  fully  extended  position.  Thus 
the  strut  does  not  begin  to  deflect  in  an  impact  until  sufficient 
force  is  developed  to  overcome  the  initial  preloading  imposed 
by  the  air  pressure  and  internal  friction.  Since  the  strut  is 
effectively  rigid  in  compression,  as  well  as  in  bending, 
prior  to  this  instant,  the  system  may  be  considered  to  have 
only  one  degree  of  freedom  during  the  initial  stage  of 


the  impact.  The  equations  of  motion  for  the  one-degree-of- 
freedom  system  are  derived  in  order  to  permit  determination 
of  the  initial  conditions  required  for  the  analysis  of  the 
landing-gear  behavior  subsequent  to  the  beginning  of  shock- 
strut  deflection. 

Since  Z\—z2—z  during  this  first  phase  of  the  impact, 
equation  (8)  maybe  written  as 

Fr,  (s)=  - y  2  ~  W(Kl-  1)  (9) 

where 


For  the  general  case  of  an  exponential  relationship  between 
vertical  ground  force  and  tire  deflection,  equation  (5)  applies 
and  the  equation  of  motion  becomes 

W 

—  z-\-mzT+W{Kh— 1)  =  0  (10) 
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The  shock  strut  begins  to  telescope  when  the  sum  of  the  inertia,  weight,  and  lift  forces  becomes  equal  to  the  vertical 
components  of  the  axial  and  normal  shock-strut  forces.  At  this  instant  tT,  Fa=Fa<)+F/r  and  equation  (6)  can  be  written  as 

(F'0+F,r)  cos  <p+F„a  sin  ip-\-KLW—  W \ 

5. - ^ -  (ID 

where 

Fa<)  initial  air-pressure  preload  force,  p„aAa 
Ffr  static  friction  at  instant  tT 

At  the  instant  tT,  s= 0  and  equation  (4)  becomes 


F,=\FNat\K> 


where 


Kf, — ^(>11  +  ^2)  + 

and  hi  and  /12  are  coefficients  of  static  friction. 

Since  the  strut  is  assumed  essentially  rigid  in  compression  (and  also  rigid  in  bending),  there  is  no  kinematic  displace¬ 
ment  of  the  lower  mass  in  the  horizontal  direction  up  to  the  beginning  of  shock-strut  deflection,  so  that  /2  =  0  and  equation 
(4b)  becomes 

/N<,t=(Vv4r+~?  zT—  W2)  sin  <f—FHgr  cos  ^  (11  b) 

Incorporating  equations  (11a),  (lib),  and  (9)  into  equation  (11)  gives 

Fa  —(±Kll  sin  v>  — cos  <p)(Kl  11’—  IT,)  —  ± 2f„cosv>+sin ip) 

— - Wt - : - 1 -  (12) 

(i-K*  sin  ip— cos  <p) 


In  equation  (12)  wherevrer  the  ±  sign  appears,  the  plus  signs 
apply  when  /V  >0  and  the  minus  signs  apply  when  F„a<0. 

From  equation  (10)  the  vertical  displacement  of  the  system 
at  the  instant  t,  is  given  in  terms  of  the  corresponding 
acceleration  by 

2'“  {^  [  w’(l *']}"  (13) 

Integrating  equation  (10)  and  noting  that  ’O=0  provides 
the  relationship  between  the  vertical  velocity  anti  the 
vertical  displacement  of  the  system  at  the  beginning  of 
shock-strut  deflection 

io2~w[}TT 2/+,+  (H) 

In  view  of  the  fact  that  the  tire  force-deflection  curve  is 
essentially  linear  for  small  deflections,  it  may  be  reasonably 
assumed  that  r=l  for  the  purpose  of  determining  the  time 
after  contact  at  which  the  strut  begins  to  telescope.  With 
this  assumption  tT  can  be  determined  from  the  relationship 


where  the  general  expression  for  the  variable  2  is  obtained 
from  equation  (14)  without  the  subscripts  r.  Performing 
the  indicated  integration  gives 

*,=V5  Ssin“  ^[(1-^)-—]}  (15) 

where 

m  S 

yj  z<?^+[(i-KL)gr 

The  computation  of  tT  can  be  greatly  simplified  by  use  of 
the  following  approximation  which  assumes  a  linear  relation¬ 
ship  between  velocity  and  time: 


1 

*r  •  1 


Equation  (15a)  should  be  a  fairly  good  approximation  in 
view  of  the  relatively  short  time  interval  between  initial  con¬ 
tact  and  the  beginning  of  shock-strut  motion. 

Equations  (12),  (13),  and  (14)  permit  the  determination 
of  the  vertical  acceleration,  displacement,  and  velocity,  re¬ 
spectively,  of  the  system  (upper  and  lower  masses)  at  the 
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beginning  of  shock-strut  deflection.  Equation  (15)  or  (15a)  permits  calculation  of  the  time  interval  between  initial  contact 
and  this  instant.  These  equations  provide  the  initial  conditions  required  for  the  analysis  of  the  behavior  of  the  landing 
gear  as  a  system  with  two  degrees  of  freedom  after  the  shock  strut  begins  to  deflect. 

If  drag  loads  are  considered,  the  solution  of  equation  (12)  requires  knowledge  of  the  horizontal  ground  force  FHs  at 
the  instant  U.  Since  the  present  analysis  does  not  explicitly  treat  the  determination  of  drag  loads,  values  of  F„s  have 
to  be  estimated,  either  from  other  analytical  considerations,  experimental  data,  or  on  the  basis  of  experience. 


MOTION  SUBSEQUENT  TO  BEGINNING  OF  SHOCK-STRUT  DEFLECTION 

Once  the  sum  of  the  inertia,  weight,  and  lift  forces  becomes  sufficiently  large  to  overcome  the  preloading  force  in  the 
shock  strut  due  to  initial  air  pressure  and  internal  friction,  the  shock  strut  can  deflect  and  the  system  becomes  one  having 
two  degrees  of  freedom.  Incorporating  the  expressions  for  the  hydraulic,  pneumatic,  and  friction  forces  (eqs.  (2a),  (3), 
and  (4))  into  equation  (6)  permits  the  equation  of  motion  for  the  upper  mass  to  be  written  as  follows: 


Wt 


Zl  + 


III 


s|  2(C, A,)2 


2+P«0a(--^"4-)  +jJ|l*Va|[(M. 


+  M,)  ^  +  *  + 


]} 


cos  Kl\V—  H’i+  Fa  sin  <^  =  0 


where 

g|  —  Sl 

COS  tp 


(16) 


and,  since  Fv’  =  FV.(^),  equation  (4c)  becomes 


COS  if 


'  v<i 


where  Fv,(z a)  is  determined  from  the  force-deflection  characteristics  of  the  tire.  For  the  ur  •  '  type  of  pneumatic  tire, 
Fv.(zt)  =  mzt,  as  previously  noted. 

Similarly,  the  equation  of  motion  for  the  lower  mass  follows  from  equation  (7): 


7s  +i|tIF-v.|[(m.+m5)^+m2]}  cos  (H-  Fr,M-Fy.  sin*-ir2=0 

The  overall  dynamic  equilibrium  equation  is  still,  of  course,  ns  given  by  equation  (8) 

V  2'  +  lI2  **+  H’(Ai-I)+*V,(a*)=0 

s  § 


(IT) 


Any  two  of  the  preceding  equations  (eqs.  (16),  (17),  and 
(8))  are  sufficient  to  describe  the  behavior  of  the  landing 
gear  subsequent  to  the  beginning  of  shock-strut  motion. 
These  equations  may  be  used  to  calculate  the  behavior  of  a 
given  landing-gear  configuration  or  to  develop  orifice  and 
metering-pin  characteristics  required  to  produce  a  specified 
behavior  for  given  impact  conditions.  They  may  also  be 
used  as  a  basis  for  the  calculation  of  dynamic  loads  in  flexible 
airplane  structures  either  by  (a)  determining  the  landing- 
gear  forcing  function  under  the  assumption  that  the  upper 
mass  is  a  rigid  body  and  then  using  this  forcing  function  to 
calculate  the  response  of  the  elastic  system  or  (b)  combining 
the  preceding  equations  with  the  equations  representing  the 
additional  degrees  of  freedom  of  the  structure;  the  simul¬ 
taneous  solution  of  the  equations  for  such  a  system  would 
then  take  into  account  the  interaction  between  the  deforma¬ 
tion  of  the  structure  and  the  landing  gear. 


SOLUTION  OF  EQUATIONS  OF  MOTION 

In  the  general  case  the  analysis  of  a  landing  gear  involves 
the  solution  of  the  equations  of  motion  given  in  the  section 
entitled  “Motion  Subsequent  to  the  Beginning  of  Shock- 
Strut  Deflection,”  with  the  initial  conditions  taken  as  the 
conditions  of  motion  at  the  beginning  of  shock-strut  deflec¬ 
tion,  as  determined  in  accordance  with  the  initial  impact 
conditions  and  the  equations  given  in  the  section  entitled 
“Motion  Prior  to  Shock-Strut  Deflection.” 

NUMERICAL  INTEGRATION  PROCEDURES 

In  view  of  the  fact  that  the  equations  of  motion  for  the 
lunding  gear  subsequent  to  the  beginning  of  shock-strut 
deflection  are  highly  nonlinear,  analytical  solution  of  these 
equations  does  not  apppar  feasible.  In  the  present  report, 
therefore,  finite-difference  methods  are  resorted  to  for  the 
step-by-step  integration  of  the  equations  of  motion.  Al- 
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though  such  numerical  methods  lack  the  generality  of  ana¬ 
lytical  solutions  and  are  especially  time  consuming  if  the 
calculations  are  carried  out  manually,  the  increasing  availa¬ 
bility  of  automatic  calculating  machines  largely  overcomes 
these  objections. 

Most  of  the  solutions  presented  in  this  report  were  obtained 
with  a  procedure,  hereinafter  referred  to  as  the  “linear  pro¬ 
cedure,”  which  assumes  changes  in  the  motion  variables  to 
be  linear  over  finite  time  intervals.  A  few  of  the  solutions 
presented  were  obtained  with  a  procedure,  hereinafter  referred 
to  as  the  “quadratic  procedure,”  which  assumes  a  quadratic 
variation  of  displacement  with  time  for  successive  intervals. 
The  generalized  solutions  for  the  simplified  equations  dis¬ 
cussed  in  a  subsequent  section  were  obtained  by  means  of 
the  Runge-Kutta  procedure.  The  application  of  these 
procedures  is  described  in  detail  in  appendix  A. 

USB  or  TIKE  FORCE-DEFLECTION  CHARACTERISTICS 

In  order  to  obtain  solutions  for  particular  cases,  it  is,  of 
course,  necessary  to  have,  in  addition  to  information  regard¬ 
ing  the  physical  characteristics  of  the  landing  gear,  some 
knowledge  of  the  force-deflection  characteristics  of  the  tire. 

If  extensive  data  regarding  the  dynamic  tire  character¬ 
istics,  such  as  shown  in  figures  2  and  3,  are  available,  an 
accurate  solution  can  be  obtained  which  takes  into  account 
the  various  breaks  in  the  force-deflection  curves  (logarithmic 
coordinates),  as  well  as  the  effects  of  hysteresis.  In  view  of 
the  fact  that  the  constants  m'  and  r  have  the  same  values 
throughout  practically  the  entire  tire  compression  process 
regardless  of  the  initial  impact  velocity  or  the  maximum 
load  attained,  these  values  of  m'  and  r,  as  determined  from 
the  force-deflection  curves,  can  be  used  in  the  calculation  of 
the  motion  subsequent  to  the  beginning  of  shock-strut  deflec¬ 
tion  until  the  first  break  in  the  force-deflection  curve  is 
reached  prior  to  the  attainment  of  the  maximum  force.  If 
the  conditions  for  the  calculations  are  the  same  as  those  for 
which  force-deflection  curves  are  available,  the  values  of 
m'  and  r  for  each  of  the  several  regimes  subsequent  to  the 
first  break  can  also  be  determined  directly  from  the  force- 
deflection  curves.  In  general,  however,  the  conditions  will 
not  be  the  same  and  interpolation  will  be  necessary  to 
estimate  the  values  of  m'  for  the  subsequent  regimes. 
Such  interpolation  is  facilitated,  particularly  after  the  maxi¬ 
mum  force-deflection  point  has  been  calculated,  by  the  fact 
that  each  subsequent  regime  has  a  fixed  value  of  r,  regardless 
of  the  initial  impact  conditions. 

The  use  of  the  tire-deflection  characteristics  in  the  calcula¬ 
tions  is  greatly  simplified  if  hysteresis  is  neglected  since  the 
values  of  m'  and  r  which  apply  prior  to  the  first  break  in  the 
force-deflection  curves  are  then  used  throughout  the  entire 
calculation,  except  in  the  case  of  severe  impacts  where  tire 
bottoming  occurs,  in  which  case  new  values  of  m ’  and  r  are 
employed  in  the  tire-bottoming  regime.  A  similar  situation 
exists  with  respect  to  the  constants  a'  and  b  when  the  linear 
approximations  which  neglect  hysteresis  are  used.  These 
simplifications  would  normally  be  employed  when  only  the 
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tire  manufacturer’s  static  or  so-called  impact  load-deflection 
data  are  available,  as  is  usually  the  case. 

EFFECT  OF  DRAG  LOADS 

Although  the  present  analysis  permits  taking  into  account 
the  effects  of  wheel  spin-up  drag  loads  on  the  behavior  of  the 
landing  gear,  the  determination  of  the  drag-load  time  history 
is  not  treated  explicitly.  Thus,  if  it  is  desired  to  consider  the 
effects  of  the  drag  load  on  the  gear  behavior,  such  as  in  the 
case  of  a  drop  test  in  which  drag  loads  are  simulated  by 
reverse  wheel  rotation  or  in  a  landing  with  forward  speed,  it 
is  necessary  to  estimate  the  drag  load,  either  by  means  of 
other  analytical  considerations  or  by  recourse  to  experimental 
data.  As  a  first  approximation  the  instantaneous  drag  force 
may  be  assumed  to  be  equal  to  the  vertical  ground  reaction 
multiplied  by  a  suitable  coefficient  of  friction  p;  that  is, 
FH  =Fv^,  up  to  the  instant  when  the  wheel  stops  skid(h>i</ 
after  which  the  drag  force  may  be  assumed  equal  to 
(The  current  ground-loads  requirements  specify  a  ski 
coefficient  of  friction  a=0.55;  limited  experimental  evidc. 
on  the  other  hand,  indicates  that  /x  may  be  as  high  as  0.7  or 
as  low  as  0.4.)  In  some  cases  experimental  data  indicate 
that  representation  of  the  drag-load  time  history  can  be 
simplified  even  further  by  assuming  a  linear  variation  of  the 
drag  force  with  time  during  the  period  of  wheel  skidding. 

The  instant  at  winch  the  wheel  stops  skidding  can  be 
estimated  from  the  simple  impulse-momentum  relationship 


Fv,  dt— 


I IT 


where 

Iu  polar  moment  of  inertia  of  wheel  assembly  about  axle 
VHq  initial  horizontal  velocity 
rd  radius  of  deflected  tire 
t,u  time  of  wheel  spin-up 

When  the  drag  force  is  expressed  in  terms  of  the  vertical 

force,  the  value  of  the  integral  j  FB  dt  can  be  determined  as 

Jo  * 

the  step-by-step  calculations  proceed  and  the  drag-force 
term  eliminated  from  the  equations  of  motion  after  the  re¬ 
quired  value  of  the  integral  at  the  instant  of  spin-up  is 
reached. 


EVALUATION  OF  ANALYSIS  BY  COMPARISON  OF  CALCULATED 
RESULTS  WITH  EXPERIMENTAL  DATA 

In  order  to  evaluate  the  applicability  of  the  foregoing 
analytical  treatment  to  actual  landing  gears,  tests  were 
conducted  in  the  Langley  impact  basin  with  a  conventional 
oleo-pneumatic  landing  gear  originally  designed  for  a  small 
military  training  airplane.  A  description  of  the  test  specimen 
and  apparatus  used  is  given  in  appendix  B. 

In  this  section  calculated  results  are  compared  with  ex¬ 
perimental  data  for  a  normal  impact  and  a  severe  impact 
with  tire  bottoming.  The  vertical  velocities  at  the  instant 
of  ground  contact  used  in  the  calculations  correspond  to 
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(he  vertical  velocities  measured  in  the  tests.  Equations 
(12),  (13),  (14),  and  (15a)  were  used  to  calculate  the  values 
of  the  variables  at  the  instant  of  initial  shock-strut  deflection. 
Numerical  integration  of  equations  (16)  and  (17)  provided 
the  calculated  results  for  the  two-degree-of-freedom  system 
subsequent  to  the  beginning  of  shock-strut  deflection. 

In  these  calculations  the  discharge  coefficient  for  the  orifice 
and  the  polytropic  exponent  for  the  air-compression  process 
were  assumed  to  have  constant  values  throughout  the  impact. 
Consideration  of  the  shape  of  the  orifice  and  examination 
of  data  for  rounded  approach  orifices  in  pipes  suggested 
a  value  of  Cd  equal  to  0.9.  Evaluation  of  data  for  other 
landing  gears  indicated  that  the  air-compression  process 
could  be  represented  fairly  well  by  use  of  an  average  value 
of  the  effective  poly  tropic  exponent  n=1.12.  In  view  of 
the  fact  that  the  landing  gear  was  mounted  in  a  vertical 
position  and  drag  loads  were  absent  in  the  tests,  friction 
forces  in  the  shock  strut  were  assumed  to  be  negligible 
in  the  calculations.  Since  the  weight  was  fully  balanced 
by  lift  forces  in  the  tests,  the  lift  factor  KL  was  taken  equal  to 
1.0.  The  appropriate  exact  tire  characteristics  (see  fig.  3) 
were  used  for  each  case. 

NORMAL  IMPACT 

Figure  4  presents  a  comparison  of  calculated  results  with 
experimental  data  for  an  impact  without  tire  bottoming  at 
a  vertical  velocity  of  8.86  feet  per  second  at  the  instant  of 
ground  contact.  The  exact  dynamic  force-deflection  charac¬ 
teristics  of  the  tire,  including  hysteresis,  were  used  in  the 
calculations.  These  tire  characteristics  are  shown  by  the 
solid  lines  in  figure  2  (a)  and  values  for  the  tire  constants 
w'  and  r  are  given  in  figure  3  (a). 

Calculated  time  histories  of  the  total  force  on  the  upper 
mass  and  the  acceleration  of  the  lower  inass  are  compared 
with  experimental  data  in  figure  4  (a).  Similar  comparisons 
for  the  upper-mass  displacement,  upper-mass  velocity,  lower- 
mass  displacement,  strut  stroke,  and  strut  telescoping 
velocity  are  presented  in  figure  4  (b).  As  can  be  seen,  the 
agreement  between  the  calculated  and  experimental  results 
is  reasonably  good  throughout  most  of  the  time  history. 
Some  of  the  minor  discrepancies  during  the  later  stages  of 
the  impact  appear  to  be  due  to  errors  in  measurement  since 
the  deviations  between  the  calculated  and  experimental 
upper-mass  accelerations  (as  represented  by  the  force  on 
the  upper  mass)  are  incompatible  with  those  for  the  upper- 
mass  displacements,  whereas  the  calculated  upper-mass  dis¬ 
placements  are  necessarily  directly  compatible  with  the 
calculated  upper-mass  accelerations.  The  maximum  value 
of  the  experimental  acceleration  of  the  lower  mass  may  be 
somewhat  high  because  of  overshoot  of  the  accelerometer. 

In  addition  to  the  total  force  on  the  upper  mass,  figure 
4  (a)  presents  calculated  time  histories  of  the  hydraulic 
and  pneumatic  components  of  the  shock-strut  force,  as 
determined  from  equations  (2)  and  (3),  respectively.  It 
can  be  seen  that  throughout  most  of  the  impact  the  force 
developed  in  the  shock  strut  arises  primarily  from  the  hy¬ 
draulic  resistance  of  the  orifice.  Toward  the  end  of  the 
impact,  however,  because  of  the  decreased  telescoping 
velocities  and  fairly  large  strokes  which  correspond  to  high 


compression  ratios,  the  air-pressure  force  becomes  larger 
than  the  hydraulic  force. 

IMPACT  WITH  TIRE  BOTTOMING 

Figure  5  presents  a  comparison  of  calculated  and  experi¬ 
mental  results  for  a  severe  impact  (VV0=  11.63  ft  per  sec) 
in  which  tire  bottoming  occurred.  The  tire  force-deflection 
characteristics  used  in  the  calculations  are  shown  by  the 
solid  lines  in  figure  3  (b).  Region  (1)  of  the  tire  force- 
deflection  curve  has  the  same  values  of  the  tire  constants 
to'  and  r  as  for  the  case  previously  discussed.  Following 
the  occurrence  of  tire  bottoming,  however,  different  values 
of  to'  and  r  apply.  These  values  are  given  in  figure  3  (b). 

It  can  be  seen  from  figure  5  that  the  agreement  between 
the  calculated  and  experimental  results  for  this  case  is 
similar  to  that  for  the  comparison  previously  presented. 

The  calculated  instant  of  tire  bottoming  is  indicated  in 
figure  5.  When  tire  bottoming  occurs,  the  greatly  increased 
stiffness  of  the  tire  causes  a  marked  increase  in  the  shock- 
strut  telescoping  velocity,  as  is  shown  in  the  right-hand 
portion  of  figure  5  (b).  Since  the  strut  is  suddenly  forced  to 
absorb  energy  at  a  much  higher  rate,  an  abrupt  increase 
in  the  hydraulic  resistance  takes  place.  The  further  increase 
in  shock-strut  force  immediately  following  the  occurrence  of 
tire  bottoming  is  evident  from  the  left-hand  portion  of 
figure  5  (a).  The  sudden  increase  in  lower-mass  acceleration 
at  the  instant  of  tire  bottoming  can  also  be  seen. 

In  this  severe  impact  the  hydraulic  resistance  of  the  orifice 
represents  an  even  greater  proportion  of  the  total  shock- 
strut  force  than  was  indicated  by  the  calculated  results  for 
an  initial  vertical  velocity  of  8.86  feet  per  second  previously 
discussed. 

The  foregoing  comparisons  indicate  that  the  analytical 
treatment  presented,  in  conjunction  with  reasonably  straight¬ 
forward  assumptions  regarding  the  parameters  involved  in 
the  equations,  provides  a  fairly  accurate  representation  of 
the  behavior  of  a  conventional  oleo-pneumatic  landing  gear. 

PARAMETER  STUDIES 

In  the  previous  section  comparisons  of  calculated  results 
with  experimental  data  showed  that  the  equations  which 
have  been  developed  provide  a  fairly  good  representation  of 
the  behavior  of  the  landing  gear  for  the  impact  conditions 
considered.  In  view  of  the  fact  that  the  equations  are 
somewhat  complicated  and  require  numerical  values  for 
several  parameters  such  as  the  tire  force-deflection  constants 
to'  and  r,  the  orifice  discharge  coefficient  Cd,  and  the  poly¬ 
tropic  exponent  «,  which  may  rot  be  readily  or  accurately 
known  in  the  case  of  practical  engineering  problems,  it 
appears  desirable  (a)  to  determine  the  relative  accuracy 
with  which  these  various  parameters  have  to  be  known  and 
(b)  to  investigate  the  extent  to  which  the  equations  can  be 
simplified  and  still  yield  useful  results.  In  order  to  accom¬ 
plish  these  objectives,  calculations  have  been  made  to 
evaluate  the  effect  of  simplifying  the  force-deflection  charac¬ 
teristics  of  the  tire,  as  well  as  to  determine  the  effects  which 
different  values  of  the  orifice  discharge  coefficient  and  the 
effective  polytropic  exponent  have  od  the  calculated  behavior. 
The  results  of  these  calculations  are  discussed  in  the  present 
section.  The  question  of  simplification  of  the  equations  of 
motion  is  considered  in  more  detail  in  a  subsequent  section. 
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(a)  Time  histories  of  forces  on  upper  mass  and  lower-mass  acceleration. 

Figure  5. — Comparisons  between  calculated  results  and  experimental  data  for  impact  with  tire  bottoming;  solution  with 
exact  exponential  tire  characteristics.  VV0  =  11.63  feet  per  second;  Cj  =  0.9;  n  =  1.12. 


(b)  Time  histories  of  landing-gear  velocities  and  displacements. 
Figure  5. — Concluded. 
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REPRESENTATION  OP  TIBS  FORCE-DEFLECTION  CHARACTERISTICS 

In  order  to  evaluate  the  degree  of  accuracy  required  for 
adequate  representation  of  the  tire  force-deflection  charac¬ 
teristics,  comparisons  are  made  of  the  calculated  behavior 
of  the  landing  gear  for  normal  impacts  and  impacts  with 
tire  bottoming  when  the  tire  characteristics  are  represented 
in  various  ways.  First,  the  force-deflection  characteristics 
will  be  assumed  to  be  exactly  as  shown  by  the  solid-line 
curves  in  figure  2  (b),  including  the  various  breaks  in  the 
curve  and  the  effects  of  hysteresis.  These  characteristics 
are  referred  to  hereinafter  as  the  exact  exponential  tire 
characteristics.  The  effects  of  simplifying  the  representa¬ 
tion  of  the  tire  characteristics  will  then  be  investigated  by 
considering  (a)  the  exponential  characteristics  without 
hysteresis;  that  is,  the  tire  will  be  assumed  to  deflect  and 
unload  along  the  same  exponential  curve,  (b)  the  linear- 
segment  approximations  to  the  tire  characteristics  (long- 
dashed  lines),  which  also  neglect  hysteresis,  and  (c)  errors 
introduced  by  neglecting  the  effects  of  tire  bottoming  in  the 
case  of  severe  impacts.  The  calculated  results  presented  in 
this  study  make  use  of  the  relationships  between  vertical 
force  on  the  tire  and  tire  deflection,  as  shown  in  figures 
3  (a)  and  3  (b). 

Figure  6  presents  a  comparison  of  the  calculated  results 
for  a  normal  impact  at  a  vertical  velocity  of  8.86  feet  per 
second,  whereas  figure  7  permits  comparison  of  the  solutions 


for  a  severe  impact,  involving  tire  bottoming,  at  a  vertical 
velocity  of  11.63  feet  per  second.  In  figures  6  and  7  the 
solid-line  curves  represent  solutions  of  the  landing-gear 
equations  when  the  exact  exponential  relationsliips  between 
force  and  tire  deflection  are  considered.  Since  these  solu¬ 
tions  were  previously  shown  to  be  in  fairly  good  agreement 
with  experimental  data  (figs.  4  and  5),  they  are  used  as  a 
basis  for  evaluating  the  results  obtained  when  tire  hysteresis 
is  neglected  and  the  force-deflection  characteristics  are  repre¬ 
sented  by  either  simplified  exponential  or  linear-segment 
relationships. 

As  in  the  calculations  previously  described,  the  solutions 
were  obtained  in  two  parts.  During  the  first  stage  of  the 
impact  the  shock  strut  was  considered  to  be  rigid  until 
sufficient  force  was  developed  to  overcome  the  initial  air- 
pressure  force.  The  calculations  for  the  landing-gear  behav¬ 
ior  subsequent  to  this  instant  were  based  on  the  equations 
which  consider  the  gear  to  have  two  degrees  of  freedom. 
Time  histories  of  the  upper-mass  acceleration  calculated  on 
the  basis  of  a  rigid  shock  strut  are  shown  by  the  dotted 
curves  in  figures  6  and  7.  These  solutions  show  the  greatest 
rate  of  increase  of  upper-mass  acceleration  possible  with 
the  exponential  tire  force-deflection  characteristics  con¬ 
sidered.  Comparison  of  these  solutions  with  those  for  the 
two-degree-of-freedom  system  indicates  the  effect  of  the 
shock  strut  in  attenuating  the  severity  of  the  impact. 


(a)  Time  histories  of  upper-mass  acceleration  and  lower-mass  acceleration. 

Figure  6. — Effect  of  tire  characteristics  on  calculated  landing-gear  behavior  in  normal  impact.  VV0=8.86  feet  per  second;  Cj=0.9;  n=1.12. 
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(a)  Time  histories  of  upper-mass  acceleration  and  lower-mass  acceleration. 

s  7. — Effect  of  tire  characteristics  on  calculated  landing-gear  behavior  for  impact  with  tire  bottoming. 
W0  =  U.fi3  feet  per  second;  Cd—O.Q;  n  =  1.12. 
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(b)  Time  histories  of  landing-gear  displacements  and  velocities. 

Figure  7.- — Continued. 
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(c)  Time  histories  of  shock-strut  stroke  and  velocity. 
Figure  7. — Concluded. 


Normal  impact. — In  the  case  of  the  normal  impact  at  a 
vertical  velocity  of  8.86  feet  per  second,  figure  6  shows  that 
the  solution  obtained  with  the  exponential  force-deflection 
variation  which  neglects  hysteresis  and  the  solution  with  the 
linear-segment  approximation  to  the  tire  characteristics  are 
in  fairly  good  agreement  with  the  results  of  the  calculation 
based  on  the  exponential  representation  of  the  exact  tire 
characteristics.  The  greatest  differences  between  the  solu¬ 
tions  are  evident  in  the  time  histories  of  upper-mass  and 
lower-mass  accelerations;  considerably  smaller  differences  are 
obtained  for  the  lower-order  derivatives,  as  might  be  ex¬ 
pected.  With  regard  to  the  upper-mass  acceleration,  the 
three  solutions  are  in  very  good  agreement  during  the  early 
stages  of  the  impact.  In  the  case  of  the  simplified  exponen¬ 
tial  characteristics,  neglect  of  the  decreased  slope  of  the 
force-deflection  curve  between  the  first  break  and  the  maxi¬ 
mum  (regime  (5)  in  fig.  3  (a))  resulted  in  the  calculation  of  a 
somewhat  higher  value  of  the  maximum  upper-mass  acceler¬ 
ation  than  was  obtained  with  the  exact  tire  characteristics. 
For  the  simplified  exponential  and  linear-segment  character¬ 
istics,  neglect  of  hysteresis  resulted  in  the  calculation  of 
somewhat  excessive  values  of  upper-mass  acceleration  sub¬ 
sequent  to  the  attainment  of  the  maximum  vertical  load. 
It  is  of  interest  to  note  that  the  calculated  results  for  the 
exponential  and  linear-segment  characteristics  without  hys¬ 
teresis  were  generally  in  quite  good  agreement  with  each 
other  throughout  the  entire  duration  of  the  impact,  although 
the  assumption  of  linear-segment  tire  force-deflection  char¬ 
acteristics  did  result  in  somewhat  excessive  values  for  the 
maximum  lower-mass  acceleration.  On  the  whole,  the 
simplified  tire  force-deflection  characteristics  considered  per¬ 
mit  calculated  results  to  be  obtained  which  represent  the 
behavior  of  the  landing  gear  in  normal  impacts  fairly  well. 


Impact  with  tire  bottoming.— In  the  case  of  the  severe 
impact  at  a  vertical  velocity  of  11.63  feet  per  second,  the 
effects  of  tire  bottoming  on  the  upper-mass  acceleration,  the 
lower-mass  acceleration,  and  the  strut  telescoping  velocity 
are  clearly  indicated  in  figure  7  by  the  calculated  results 
based  on  the  exact  tire  characteristics.  As  can  be  seen,  the 
linear-segment  approximation  to  the  tire  deflection  character¬ 
istics  which  takes  into  account  the  effects  of  tire  bottoming 
resulted  in  a  reasonably  good  representation  of  the  landing- 
gear  behavior  throughout  most  of  the  time  history.  On  the 
other  hand,  as  might  be  expected,  the  calculations  winch 
neglected  the  effects  of  bottoming  on  the  tire  force-deflection 
characteristics  did  not  reveal  the  marked  increase  in  the 
upper-mass  acceleration  due  to  the  increased  stiffness  of  the 
tire  subsequent  to  the  occurrence  of  bottoming.  It  is  also 
noted  that  the  discrepancies  in  the  calculated  upper-mass 
acceleration  due  to  neglect  of  hysteresis  in  the  later  stages 
of  the  impact  are  more  pronounced  in  this  case  than  in  the 
impact  without  tire  bottoming  previously  considered,  again 
as  might  be  expected. 

effect  of  orifice  discharge  coefficient 

In  view  of  the  fact  that  there  is  very  little  information 
available  regarding  the  magnitude  of  discharge  coefficients 
for  orifices  in  landing  gears,  it  appears  desirable  to  evaluate 
the  effect  which  differences  in  the  magnitude  of  the  orifice 
coefficient  can  have  on  the  calculated  results.  Figure  8 
presents  comparisons  of  calculated  results  for  a  range  of 
values  of  the  orifice  discharge  coefficient  Cd  between  1.0 
and  0.7.  The  four  solutions  presented  are  for  the  same  set 
of  initial  conditions  as  the  normal  impact  without  tire 
bottoming  previously  considered  and  are  based  on  the 
exponential  tire  force-deflection  characteristics  which  neglect 
hysteresis. 
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(a)  Time  histories  of  upper-mass  acceleration  and  lower-mass  acceleration. 

Figure  8. — Effect  of  orifice  discharge  coefficient  on  landing-gear  behavior;  calculations  with  exponential  tire  characteristics 

without  hysteresis.  VV0=8.86  feet  per  second;  n=1.12. 


(b)  Time  histories  of  landing-gear  displacements  and  velocities. 
Figure  8. — Continued. 


287846—5. 


22 


REPORT  1154 — NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


(c)  Time  histories  of  shock-strut  stroke  and  velocity. 
Figure  8. — Concluded. 


Those  calculations  show  that  a  decrease  in  the  orifice  dis¬ 
charge  coefficient  results  in  an  approximately  pioportional 
increase  in  the  upper-mass  acceleration.  This  vari¬ 
ation  is  to  he  expected  since  the  smaller  coefficients  cor¬ 
respond  to  reduced  effective  orifice  areas  which  result  in 
greater  shock-strut  forces  due  to  increased  hydraulic 
resistance.  As  a  result  of  the  increased  shock-strut  force 
acting  downward  on  the  lower  mass,  the  maximum  upward 
acceleration  of  the  lower  mass  is  reduced  with  decreasing 
values  of  the  discharge  coefficient.  The  increase  in  shock- 
strut  force  with  decreasing  discharge  coefficient  also  results 
in  a  decrease  in  the  strut  stroke  and  telescoping  velocity  but 
an  increase  in  the  lower-mass  velocity  and  displacement, 
as  might  be  expected.  However,  since  the  increases  in 
lower-mass  displacement  and  velocity  are  smaller  than  the 
decreases  in  strut  stroke  and  telescoping  velocity,  the  upper- 
mass  displacement  and  velocity  are  reduced  with  decreasing 
orifice  discharge  coefficient. 

These  comparisons  show  that  the  magnitude  of  the  orifice 
coefficient  has  an  important  effect  on  the  behavior  of  the 
landing  gear  and  indicates  that  a  fairly  accurate  determi¬ 
nation  of  the  numerical  value  of  this  parameter  is  necessary 
to  obtain  good  results. 

EFFECT  OF  AU-COMPBESSION  PROCESS 

Since  the  nature  of  the  air-compression  process  in  a  shock 
strut  is  not  well-defined  and  different  investigators  have 
assumed  values  for  the  poly  tropic  exponent  ranging  any¬ 
where  between  the  extremes  of  1.4  (adiabatic)  and  1.0 


(isothermal),  it  appeared  desirable  to  evaluate  the  im¬ 
portance  of  the  air-compression  process  and  to  determine  the 
extent  to  which  different  values  of  the  polytropic  exponent 
can  influence  the  calculated  results.  Consequently,  solutions 
have  been  obtained  for  three  different  values  of  the  poly- 
tropic  exponent;  namely,  n  =  1.3,  1.12,  and  0. 

The  value  «  =  1.3  corresponds  to  a  very  rapid  compression 
in  which  an  adiabatic  process  is  almost  attained.  The 
value  n  =  1.12  corresponds  to  a  relatively  slow  compression 
in  which  the  process  is  virtually  isothermal.  The  value 
»!=0  is  completely  fictitious  since  it  implies  constant  air 
pressure  within  the  strut  throughout  the  impact.  The 
assumption  v  —  0  has  been  considered  since  it  makes  one  of 
the  terms  in  the  equations  of  motion  a  constant  and  permits 
simplification  of  the  calculations.  The  three  solutions 
presented  are  for  the  same  set  of  initial  conditions  as  the 
normal  impact  without  tire  bottoming  previously  con¬ 
sidered  and  are  based  on  the  exponential  tire  force-deflection 
characteristics  which  neglect  hysteresis. 

Figure  9  shows  that  the  air  pressure  contributes  only  a 
relatively  small  portion  of  the  total  shock-strut  force  through¬ 
out  most  of  the  impact  since  the  compression  ratio  is  rela¬ 
tively  small  until  the  later  stages  of  the  impact.  Toward 
the  end  of  the  impact,  however,  the  air-pressure  force 
becomes  a  large  part  of  the  total  force  since  the  compression 
ratio  becomes  large,  whereas  the  hydraulic  resistance  de¬ 
creases  rapidly  as  the  strut  telescoping  velocity  is  reduced 
to  zero. 
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(c)  Time  histories  of  shock-strut  stroke  and  velocity. 
Figure  9.- — Concluded. 


As  a  result,  the  calculations  show  that  the  magnitude  of 
the  poly  tropic  exponent  has  only  a  very  small  effect  on  the 
behavior  of  the  landing  gear  throughout  most  of  the  impact. 
For  the  practical  range  of  polytropic  exponents,  variations 
in  the  air-compression  process  result  in  only  minor  differ¬ 
ences  in  landing-gear  behavior,  even  during  the  very  latest 
stages  of  the  impact.  The  assumption  of  constant  air 
pressure  in  the  strut  throughout  the  impact  (n=0),  however, 
does  lead  to  the  calculation  of  excessive  values  of  stroke  and 
of  the  time  to  reach  the  maximum  stroke.  The  time  history 
of  the  shock-strut  force  calculated  on  the  basis  of  this 
assumption  is,  on  the  other  hand,  in  quite  good  agreement 
with  the  results  for  the  practical  range  of  air-compression 
processes. 

On  the  whole  it  appears  that  the  behavior  of  the  landing 
gear  is  relatively  insensitive  to  variations  in  the  air- 
compression  process.  The  foregoing  results  suggest  that,  in 
many  cases,  fairly  reasonable  approximations  for  the  landing- 
gear  force-time  variation  might  be  obtained  even  if  the  air- 
pressure  term  in  the  equations  of  motion  were  completely 
neglected. 

SIMPLIFICATION  OF  EQUATIONS  OF  MOTION 

The  preceding  studies  have  indicated  that  variations  in  the 
tire  force-deflection  characteristics  and  in  the  air-compression 
process  individually  have  only  a  relatively  minor  effect  on  the 
calculated  behavior  of  the  landing  gear.  These  results  sug¬ 
gest  that  the  equations  of  motion  for  the  landing  gear  might 
be  simplified  by  completely  neglecting  the  internal  air- 
pressure  forces  in  the  shock  strut  and  by  considering  the  tiro 
force-deflection  characteristics  to  be  linear.  With  these  as¬ 
sumptions,  the  equations  of  motion  for  the  upper  mass, 
lower  mass,  and  complete  system  (eqs.  (16),  (17),  and  (8))  can 


be  written  as  follows  for  the  case  where  the  wing  lift  is  equal 
to  the  weight  and  the  internal  friction  is  neglected: 


+  — i2)J-)-WT2=0 


Zz — A(ii — Z2)24-ttS2'4~ b — \I2~0  ► 

Wt  ..  ,Wt.  ,  . 

- 2 1  — 1 - 02~b  IZ2A  0  =  0 

g  g 

where 


(18) 


A=i-  pAh3 
|s|  2(C„An)2cosifi 

and 

a  slope  of  linear  approximation  to  tire  force-deflection 
characteristics 

b  value  of  force  corresponding  to  zero  tire  deflection,  as 
determined  from  the  linear-segment  approximation  to 
the  tire  force-deflection  characteristics 
The  motion  variables  at  the  beginning  of  shock-strut  de¬ 
flection  can  be  readily  determined  in  a  manner  similar  to 
that  employed  in  the  more  general  treatment  previously  dis¬ 
cussed.  For  the  simplified  equations  the  variables  at  the 
instant  tT  are  given  by 


Z  T= 


Wj 

Wx 


g 


(19) 


Sr 


na 

In  most  cases  the  term  Zr2  is  small  in  comparison  with  i02 
so  that  iT«z0. 
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(a)  Time  histories  of  upper-mass  acceleration  and  lower-mass  acceleration. 
Figure  10. — Evaluation  of  calculated  results  for  simplified  systems.  Vr0=8.86  feet  per  second;  C*— 0.9. 


The  values  determined  from  equations  (19)  are  used  as 
initial  conditions  in  the  solution  of  equations  (18). 

The  fact  that  the  lower  mass  is  a  relatively  small  fraction  of 
the  total  mass  suggests  that  the  system  might  be  simplified 
even  further  without  greatly  modifying  the  calculated  results 
by  assuming  the  lower  mass  to  be  equal  to  zero.  With  this 
assumption  f,=0  and  the  initial  values  of  the  variables  in 
equations  (18)  correspond  to  the  conditions  at  initial  contact. 

EVALUATION  OF  SIMPLIFICATION 

In  order  to  evaluate  the  applicability  of  these  simplifica¬ 
tions,  the  behavior  of  the  landing  gear  has  been  calculated  in 
accordance  with  equations  (18)  for  an  impact  with  an  initial 
vertical  velocity  of  8.86  feet  per  second.  A  similar  calcula¬ 
tion  has  been  made  with  the  assumption  lfr2=0.  These 
results  are  compared  in  figure  10  with  the  more  exact  solu¬ 
tions  previously  presented  in  figure  4,  which  include  con¬ 
sideration  of  the  air-compression  springing  and  the  exact 
exponential  tire  characteristics.  A  time  history  of  the  lower- 
mass  acceleration  is  not  presented  for  the  case  where  W 2  is 
assumed  equal  to  zero  since  the  values  of  z2jg  have  no 
significance  in  this  case. 

Figure  10  shows  that  the  two  simplified  solutions  are  in 
quite  good  agreement  with  each  other,  as  might  be  expected, 
and  are  also  in  fairly  good  agreement  with  the  more  exact 
results.  Neglecting  the  air-pressure  forces  and  assuming  a 
linear  tire  force-deflection  variation  resulted  in  the  calcula¬ 
tion  of  slightly  lower  values  for  the  maximum  upper-mass 
acceleration  and  somewhat  higher  values  for  the  maximum 
stroke  than  were  obtained  with  the  more  exact  equations. 


The  effect  of  neglecting  the  lower  mass  was  primarily  to 
reduce  the  lower-mass  displacement  (tire  deflection),  as  a 
result  of  the  elimination  cr  the  lower-mass  inertia  reaction. 

On  the  whole,  it  appea  that  the  assumptions  considered 
permit  appreciable  simpli'"  nation  of  the  equations  of  motion 
without  greatly  impairing  *  :  validity  of  the  calculated  results. 

GENERALIZED  TREATMENT 

Equations  and  solutions. — By  writing  the  simplified  equa¬ 
tions  of  motion  in  terms  of  dimensionless  variables,  general¬ 
ized  solutions  can  be  obtained  for  a  wide  range  of  landing-gear 
and  impact  parameters  which  may  be  useful  in  pre¬ 
liminary  design.  If  W2  is  taken  equal  to  zero  and  it  is 
further  assumed  that  the  tire  force-deflection  curve  is 
represented  by  a  single  straight  line  through  the  origin 
(6=0  throughout  the  impact),  equations  (18)  reduce  to 

'±it+An,-i,r=o ) 

A(zi— z2)2—az2=0  >  (20) 

Wi.  ,  n  \ 

—  zi+az2=0  I 

g  / 

W, 

where  — ->  A,  and  a  are  constants,  as  previously  defined,  and 

any  two  of  the  foregoing  equations  are  sufficient  to  describe 
completely  the  behavior  of  the  system .  With  this  representa¬ 

tion  of  the  system,  the  shock  strut  begins  to  deflect  at  the  in¬ 
stant  of  initial  contact  (f,=0).  Thus,  the  initial  conditions 
for  equations  (20)  are  the  initial  impact  conditions;  namely, 
2>0=2j0=0  and  2I0=i2o=2°- 
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As  can  be  seen,  with  the  equations  in  this  form,  the  solution 

W 

depends  on  five  parameters,  namely,  A,  a,  and  the 

initial  conditions  i0  and  So.  However,  since  z0=0  in  all 
cases,  the  number  of  variable  parameters  is  reduced  to  four. 
In  view  of  the  fact  that  these  parameters  are  independent  of 
one  another  and  each  may  take  on  a  large  range  of  values,  a 
great  many  solutions  and  a  large  number  of  graphs  would  be 
required  to  cover  the  entire  range  of  landing-gear  and  impact 
parameters  with  the  equations  in  the  form  of  equations  (20). 

The  number  of  independent  parameters  which  have  to  be 
considered  may  be  greatly  decreased  by  the  introduction  of 
generalized  dimensionless  variables  and  the  corresponding 
transformations  of  equations  (20).  In  this  case,  generalized 
variables  can  be  obtained  which  permit  transformation  of  the 
equations  of  motion  to  a  form  which  does  not  involve  any 
constants.  With  the  equations  in  this  form,  there  is  only 
one  variable  parameter,  namely,  the  initial  velocity  param¬ 
eter.  To  determine  the  generalized  variables  which  satisfy 
the  aforementioned  requirements,  let 


and 

0=/ 

Thus, 

,  du 

U  ~dP 

and 

du' 

u  ~do 

Substituting  these  new  variables  permits  equations  (20) 
to  be  written  as 


(it/-it/)2+(-|  u,"=0 

(it/-it/)2— «2=0 

u"+(ma) U2=0 


(20a) 


The  number  of  independent  parameters  will  be  reduced 
if  all  the  combined  constants  in  equations  (20a)  are  set 
equal  to  one  another,  that  is,  let 


Wt _ a  a  _  alp2 

A  g-A?-WTt 


From  this  relationship,  it  can  be  seen  that 

A 


and 


Wt/g 


0 


~ylwxig 


Thus  the  generalized  variables  become 


~Zl(wt/g) 


„  ,  du ,  .  I  A2 1  a 

Ul -d8==z'yJwTg 


Ut~Zt{w^) 


U2~de-Z2 


A2/  a 
WJg 


and 

„  d2ui  .  ( A\ 
u'-dei=z'Ka) 

where 

1  a 

With  these  new  variables  equations  (20) 

(it/  —  Mj')*+«i"  =  0\ 

(it/  — Itj')*— «2=0  > 

it/'+ita=0  ) 

(21) 


where  any  two  of  these  equations  are  sufficient  to  describe 
the  behavior  of  the  system. 

Inasmuch  as  equations  (21)  do  not  involve  any  constants, 
their  solutions  are  completely  determined  by  the  initial  val¬ 
ues  of  the  variables.  Since  the  displacements  at  initial 
contact  ut()  and  u2q  are  equal  to  zero  and  the  initial  velocities 
iti0'  and  u^'  are  equal,  the  only  parameter  is  the  initial 
dimensionless  velocity 


It/  —  Zq 


V  Wxa 


where  it/ = it,/ =112/- 

Generalized  solutions  of  equations  (21)  are  presented  in 
figure  11  for  values  of  it/  corresponding  to  a  wide  range  of 
landing-gear  and  impact  parameters.  Parts  (a)  to  (e)  of 
figure  11  show  the  variations  of  the  dimensionless  variables 
during  the  impact;  parts  (f)  and  (g)  show  the  maximum 
values  of  the  more  important  variables  as  functions  of  it/. 
Part  (h)  shows  the  shock-strut  effectiveness  ij,  and  the 
landing-gear  effectiveness  The  shock-strut  effective¬ 

ness,  sometimes  called  “efficiency”  and,  in  Europe,  “plani- 
metric  ratio,”  is  defined  as 

J'Vma, 

uTda 


where  o-=m1— u2  is  the  dimensionless  shock-strut  stroke. 
Since  ij,  represents  the  ratio  of  the  energy  actually  absorbed 
by  the  shock  strut  to  the  maximum  energy  which  the  strut 
could  possibly  absorb  for  any  combination  of  maximum 
acceleration  (or  load)  and  maximum  stroke,  it  serves  as  a 
measure  of  the  extent  to  which  a  given  combination  of 
maximum  load  and  stroke  has  been  utilized  to  absorb  the 
energy  of  an  impact.  A  similar  measure  of  the  energy 
absorption  effectiveness  of  the  landing  gear  as  a  whole  is 
given  by  rn,  which  is  defined  by 


nt.= 


J’Uj 


ti,"  dut 


i*i 


*  Equations  (21)  may  be  reduced  to  a  single  equation  in  one  variable  by  differentiating  the 
last  equation  and  substituting  for  tu'  in  the  first  equation.  This  gives 
By  introducing  the  new  variable  this  equation  may  be  reduced  to  the  second- 

order  equation  (te+*c"),+w'”0,  subject  to  the  initial  conditions  and  *V*0. 
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(a)  Relationship  between  upper-mass  acceleration,  lower-mass  displacement,  and  time. 

Figure  11. — Generalized  solutions  for  simplified  system. 


(b)  Relationship  between  upper-mass  displacement  and  time. 
Figure  11. — Continued. 
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(c)  Relationship  between  upper-mass  velocity,  lower-mass  velocity,  and  time. 
Figure  11. — Continued. 


The  generalized  results  presented  in  figure  11 3  can  be 
used  to  estimate  the  performance  of  a  given  landing  gear  of 
known  configuration  for  particular  impact  conditions  or  to 
choose  the  dimensions  for  a  landing  gear  when  the  impact 
conditions  and  desired  performance  are  specified. 

Applicability  of  solutions. — To  illustrate  the  applicability 
of  the  generalized  solutions,  the  curves  of  figure  11  have 
been  applied  to  the  previously  considered  case  of  the  normal 
impact  at  an  initial  vertical  velocity  of  8.86  feet  per  second 
for  comparison  with  the  more  exact  solution  presented  in 
figure  4.  In  order  to  make  use  of  the  generalized  solutions 

*  Although  time-history  solutions  are  presented  for  values  of  ««'  as  small  as  0.5,  It  will  be 
noted  that  values  of  ulmS0t  9mm,  ?»,  and  s/r  are  not  given  for  values  of  u«'<1.5.  It  can  be 
seen  from  the  time  histories  tl»at  the  characteristics  of  the  solution  in  the  later  stages  of  the  im¬ 
pact  change  as  u%'  becomes  small;  in  particular,  ut  increases  and  the  curve  of  ftinM  as  a  func¬ 
tion  of  W  appears  to  reach  a  minimum  at  some  value  of  u«'31.A.  Furthermore,  the  later 
stages  of  the  solutions  greatly  stretch  out  in  time  and  appear  to  be  almost  asymptotic  in 
character.  Several  different  analytical,  numerical,  and  analogue  methods  were  applied  in  an 
attempt  to  study  this  phase  of  the  problem  further  but  the  extremely  slow  rate  of  change  of 
the  variables  in  this  region  prevented  successful  completion  of  the  solutions. 

287846 — 54 - 4 


it  is  first  necessary*  to  approximate  the  tire  force-deflection 
characteristics  by  a  simple  linear  variation.  Two  such 
linear  approximations  which  might  be  considered  suitable 
for  this  purpose  are  shown  in  figure  12.  Linear  approxi¬ 
mation  I  is  a  straight  line  through  the  origin  having  a  slope 
a=18.5X103  pounds  per  foot  (a'=ad=41.6X  103  lb).  This 
value  of  a  and  the  other  pertinent  landing-gear  and  impact 
parameters  result  in  a  value  of  the  initial  dimensionless 
velocity  parameter  u0,=2.57.  Linear  approximation  II  is 
a  straight  line  with  slope  a=21.3X103  pounds  per  foot 
(a'=47.9X  103  lb)  which  does  not  pass  through  the  origin 
but  intersects  the  displacement  axis  at  a  value  of 
Zj  =0.0508  foot.  With  this  value  of  a,  tt/— 2.39. 

(  vi“V 

Since  the  solutions  of  figure  11  have  been  calculated  only 
for  integral  values  of  Uo',  curves  for  the  foregoing  values  of 
«o'  were  graphically  interpolated  by  cross  plotting.  These 
results  were  then  converted  to  dimensional  values  by  multi- 
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(d)  Relationship  between  upper-mass  acceleration  and  strut  stroke. 
Fiouhe  11. — Continued. 


plying  the  dimensionless  variables  by  the  appropriate 
constants.  The  results  obtained  are  compared  in  figure  13 
with  the  more  exact  solution  presented  in  figure  4.  The 
values  based  on  linear  approximation  I  have  been  plotted 
exactly  as  determined  from  the  generalized  solutions.  The 
results  for  linear  approximation  II,  however,  have  been 
displaced  relative  to  the  origin  of  coordinates  as  indicated 
in  the  following  discussion. 

The  assumption  of  linear  approximation  II  implies  that 


the  system  must  move  a  distance  equal  to  z2 


after 


initial  contact  (at  constant  velocity  since  the  wing  lift  is 
taken  equal  to  the  weight)  before  any  finite  ground  reaction 
can  develop.  The  derivation  of  the  equations  of  motion, 
on  the  other  hand,  assumes  that  the  ground  reaction  in¬ 
creases  linearly  with  deflection  from  the  instant  of  initial 
contact.  As  a  result,  the  equations  of  motion  do  not  apply 
until  after  the  system  has  attained  a  displacement  equal  to 
z2(r  _oV  which  occurs  at  a  time  after  initial  contact 

(  vt~  ) 


2j 


W,-°) 


In  other  words,  the  equations  of  motic 

■  r« 

apply  to  a  coordinate  system  transformed  so  that  the  til 

force-deflection  relationship  passes  through  the  origin;  thi 

is,  a  coordinate  system  displaced  by  z%  relative  t 

(  vt~v) 


the  coordinate  system  originating  at  the  point  of  initial 
contact.  It  therefore  follows  that  the  upper-mass  and 
lower-mass  displacements  determined  from  the  generalized 
solutions  for  the  case  of  linear  approximation  II  must  be 


increased  by  a  constant  amount  equal  to  zt. 


ease  0.0508  foot,  and  all  results  must  be  displaced 


in  this 


in 


~2(zv  -o) 

time  bv  a  constant  increment  A t=  \T  —  >  in  this  case 


0.0508 

8.86 


Vr„ 

=0.0057  second,  relative  to  the.  instant  of  initial 


contact.  These  corrections  have  been  incorporated  in  plot¬ 
ting  the  curves  for  linear  approximation  II  shown  in  figure  13. 

As  can  be  seen,  the  results  obtained  by  application  of  the 
generalized  solutions,  particularly  by  the  method  employing 
linear  approximation  II,  are  in  fairly  good  agreement  with  the 
more  exact  solution.  The  discrepancies  which  exist  are 
attributable  to  the  neglect  of  the  shock-strut  preloading  and 
springing  provided  by  the  air-pressure  force,  neglect  of  the 
lower  mass,  and  to  differences  between  the  very  simple  tire 
force-deflection  relationships  assumed  and  the  exact  tire 
characteristics.  On  the  whole,  it  appears  that  the  general¬ 
ized  results  offer  a  means  for  rapidly  estimating  the  behavior 
of  the  landing  gear  within  reasonable  limits  of  accuracy  and 
may  therefore  be  useful  for  preliminary  design  purposes. 
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Omcnaonttii  uppar-mou  ditploctmant,  </, 

(e)  Relationship  between  upper-mass  acceleration  and  upper-mass  displacement. 
Fic.urk  11. — Continued. 


SUMMARY  OF  RESULTS  AND  CONCLUSIONS 

A  theoretical  study  has  been  made  of  the  behavior  of  the 
conventional  type  of  oleo-pneumatic  landing  gear  during  the 
process  of  landing  impact.  The  basic  analysis  is  presented 
in  a  general  form  and  treats  the  motions  of  the  landing  gear 
prior  to  and  subsequent  to  the  beginning  of  shock-strut 
deflection.  In  the  first  phase  of  the  impact  the  landing  gear 
is  treated  as  a  single-degree-of-freedom  system  in  order  to 
determine  the  conditions  of  motion  at  the  instant  of  initial 
shock-strut  deflection,  after  which  instant  the  landing  gear  is 
considered  as  a  system  with  two  degrees  of  freedom.  The 
equations  for  the  two-degree-of-freedom  system  consider 
such  factors  as  the  hydraulic  (velocity  square)  resistance  of 
the  orifice,  the  forces  due  to  air  compression  and  internal 
friction  in  the  shock  strut,  the  nonlinear  force-deflection 
characteristics  of  the  tire,  the  wing  lift,  the  inclination  of  the 
landing  gear,  and  the  effects  of  wheel  spin-up  drag  loads. 

The  applicability  of  the  analysis  to  actual  landing  gears 
has  been  investigated  for  the  particular  case  of  a  vertical 
landing  gear  in  the  absence  of  drag  loads  by  comparing 
calculated  results  with  experimental  drop-test  data  for  corre¬ 
sponding  impact  conditions,  for  both  a  normal  impact  and  a 
severe  impact  involving  tire  bottoming. 

Studies  have  also  been  made  to  determine  the  effects  of 
variations  in  such  parameters  as  the  dynamic  force-deflection 


characteristics  of  the  tire,  the  orifice  discharge  coefficient,  and 
the  effective  poly  tropic  exponent  for  the  air-compression 
process,  which  might  not  be  known  accurately  in  practical 
design  problems. 

In  addition  to  the  more  exact  treatment  an  investigation 
has  also  been  made  to  determine  the  extent  to  which  the 
basic  equations  of  motion  can  be  simplified  and  still  yield 
useful  results.  Generalized  solutions  of  the  simplified 
equations  obtained  are  presented  for  a  wide  range  of  landing- 
gear  and  impact  parameters. 

On  the  basis  of  the  foregoing  studies  the  following  con¬ 
clusions  are  indicated: 

1.  The  behavior  of  the  landing  gear  as  calculated  from  the 
basic  equations  of  motion  was  found  to  be  in  good  agreement 
with  experimental  drop-test  data  for  the  case  of  a  vertical 
landing  gear  in  the  absence  of  drag  loads,  for  both  a  normal 
impact  and  a  severe  impact  involving  tire  bottoming. 

2.  A  study  of  the  effects  of  variations  in  the  force-deflection 
characteristics  of  the  tire  indicates  that 

a.  In  the  case  of  a  normal  impact  without  tire  bottoming, 
reasonable  variations  in  the  force-deflection  characteristics 
of  the  tire  have  only  a  relatively  small  effect  on  the  calculated 
behavior  of  the  landing  gear.  Approximating  the  rather 
complicated  force-deflection  characteristics  of  the  actual  tire 
by  simplified  exponential  or  linear-segment  variations  appears 
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(f)  Variation  of  maximum  upper-mass  acceleration  and  time  to  reach  maximum  upper-mass  acceleration  with  initial  velocity  parameter. 

Figure  11. — Continued. 


to  be  adequate  for  practical  purposes.  Tire  hysteresis  was 
found  to  be  relatively  unimportant. 

b.  In  the  case  of  a  severe  impact  involving  tire  bottoming, 
the  use  of  simplified  exponential  and  linear-segment  approxi¬ 
mations  to  the  actual  tire  force-deflection  characteristics 
which  neglect  the  effects  of  tire  bottoming,  although  adequate 
up  to  the  instant  of  bottoming,  fails  to  indicate  the  pro¬ 
nounced  increase  in  landing-gear  load  which  results  from 
bottoming  of  the  tire.  The  use  of  exponential  or  linear- 
segment  approximations  to  the  tire  characteristics  which 
take  into  account  the  increased  stiffness  of  the  tire  that  re¬ 
sults  from  bottoming,  however,  yields  good  results. 

3.  A  study  of  the  importance  of  the  discharge  coefficient 
of  the  orifice  indicates  that  the  magnitude  of  the  discharge 
coefficient  has  a  marked  effect  on  the  calculated  behavior  of 
the  landing  gear;  a  decrease  in  the  discharge  coefficient  (or 
the  product  of  the  discharge  coefficient  and  the  net  orifice 
area)  results  in  an  approximately  proportional  increase  in 
the  maximum  upper-mass  acceleration. 

4.  A  study  of  the  importance  of  the  air-compression  process 
in  the  shock  strut  indicates  that  the  air  springing  is  of  only 
minor  significance  throughout  most  of  the  impact,  and  that 
variations  in  the  effective  polytropic  exponent  n  between  the 
isothermal  value  of  1.0  and  the  near-adiabatic  value  of  1.3 
have  only  a  secondary  effect  on  the  calculated  behavior  of 


the  landing  gear.  Even  the  assumption  of  eonstant  air 
pressure  in  the  strut  equal  to  the  initial  pressure  (n=0) 
yields  fairly  good  results,  which  may  be  adequate  for  many 
practical  purposes. 

5.  An  investigation  of  the  extent  to  which  the  equations  of 
motion  for  the  landing  gear  can  be  simplified  and  still  yield 
acceptable  calculated  results  indicates  that,  for  many  prac¬ 
tical  purposes,  the  air-pressure  force  in  the  shock  strut  can 
be  completely  neglected,  the  tire  force-deflection  relationship 
can  be  assumed  to  be  linear,  and  the  lower  or  unsprung  mass 
can  be  taken  equal  to  zero. 

6.  Generalization  of  the  equations  of  motion  for  the 
simplified  system  described  in  the  preceding  paragraph 
shows  that  the  behavior  of  this  system  is  completely  deter¬ 
mined  by  the  magnitude  of  one  parameter,  namely,  the 
dimensionless  initial-velocity  parameter.  Solution  of  these 
generalized  equations  in  terms  of  dimensionless  variables 
permits  compact  representation  of  the  behavior  of  the  system 
for  a  wide  range  of  landing-gear  and  impact  parameters, 
which  may  be  useful  for  rapidly  estimating  landing-gear 
performance  in  preliminary  design. 

Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  May  1, 1962. 
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Linear  approximation  I:-„  /  /  /  Experimental: 

<7=I85XI03  Ib/ft  'y  i  /  VVo  >  8.86  fps 


/  . 


0 


/  \ 


//L 

/  /"  Linear  approximation  I: 
/  <7  =  2I.3XI03  Ib/ft 


- Solution  of  figure  4 

Generalized  solutions 

- u0'  ■  2.57  (<7  *  I8.5XI03) 

- u‘  -  2.39  (O-2I.3XI03) 


Tire  deflection,  ft 


M  h-A/  .04  08  .12  .16  .20 

Time  ofter  contact,  sec 


Figure  12. — Linear  approximations  to  tire  force-deflection  character-  (a)  Time  history  of  upper-mass  acceleration. 

istics  used  in  application  of  generalized  solutions.  Figure  13,-Comparison  of  generalized  results  and  more  exact  solu 

tion.  1  r0  =  8.8C  feet  per  second:  Cj  —  0. 9. 


Strut  stroke,  ft 


APPENDIX  A 

NUMERICAL  INTEGRATION  PROCEDURES 


As  previously  noted,  most  of  the  specific  solutions  presented 
in  this  report  were  obtained  with  a  numerical  integration 
procedure,  termed  the  “linear  procedure,”  which  assumes 
changes  in  the  variables  to  be  linear  over  finite  time  intervals. 
With  this  procedure  a  time  interval  *  =  0.001  second  was 
used  in  order  to  obtain  the  desired  accuracy  for  the  particular 
cases  considered.  A  few  of  the  specific  solutions  presented 
were  obtained  by  means  of  a  procedure,  termed  the  “quad¬ 
ratic  procedure,”  which  assumes  a  quadratic  variation  of 
displacement  with  time  for  successive  intervals.  This  pro¬ 
cedure,  although  requiring  somewhat  more  computing  time 
per  interval,  may  permit  an  increase  in  the  interval  size  for 
a  given  accuracy  in  some  cases  allowing  a  reduction  in  the 
total  computing  time  required.  In  the  case  of  the  more 
exact  equations  of  motion  the  accuracy  of  the  quadratic  pro¬ 
cedure  with  a  time  interval  of  0.002  second  appears  to  be 
equal  to  that  of  the  linear  procedure  with  an  interval  of 
0.001  second.  Although  the  accuracy  naturally  decreases 
with  increasing  interval  size,  the  loss  in  accuracy  for  pro¬ 
portionate  increases  in  interval  size  appears  to  lie  smaller  for 
the  quadratic  than  for  the  linear  procedure.  In  the  case  of 
the  simplified  equations  of  motion  reasonably  satisfactory 
results  were  obtained  in  test  computations  with  the  quadratic 
procedure  for  intervals  as  large  as  0.01  second,  whereas  the 
linear  procedure  was  considered  questionable  for  intervals 
larger  than  0.002  second. 

The  generalized  solutions  presented,  because  of  the 
relatively  simple  form  of  the  equations  of  motion,  were 
obtained  with  the  well-known  Runge-Kutta  procedure, 
A  study  of  the  allowable  interval  size  resulted  in  the  use  of 
an  interval  A0  =  O. 08,  which  corresponds  to  a  time  interval  of 
about  0.00.5  second  for  the  landing  gear  under  consideration. 

LINEAR  PROCEDURE 

In  this  step-by-step  procedure  the  variations  in  dis¬ 
placement,  velocity,  and  acceleration  are  assumed  to  be 
linear  over  each  finite  time  interval  *.  The  metiiod,  as 
used,  involves  one  stage  of  iteration.  Linear  extrapolation 
of  the  velocity  at  the  end  of  any  interval  is  used  to  obtain 
estimated  values  of  velocity  and  displacement  for  the  next 
interval.  These  values  are  then  used  to  calculate  values  of 
the  acceleration  in  accordance  with  the  equations  of  motion. 
Integration  of  the  acceleration  provides  improved  values  of 
the  velocity  and,  if  desired,  the  displacement  and  accelera¬ 
tion.  In  this  procedure  all  integrations  are  performed  by 
application  of  the  trapezoidal  rule. 

The  following  derivation  illustrates  the  application  of  the 
linear  procedure  to  the  equations  of  motion  for  the  landing 
gear,  which  apply  subsequent  to  the  beginning  of  shock-strut 
deflection  at  time  t„  In  the  example  presented  internal 
friction  forces  are  neglected  in  order  to  simplify  the  deriva- 
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lion.  However,  the  same  general  procedure  can  be  used  if 
these,  or  other  complicating  effects,  arc  included  in  the 
equations. 

For  the  case  under  consideration  the  equations  of  motion 
(eqs.  (16),  (17).  and  (8))  can  be  written  as  follows: 

y-  21+A(i1-i2)4+«[l-t,(;,-.-2)l-”  =  />  =  0  (A  1 ) 

~  St-Afr-ztY-mi-CiZt-ztn-'  +  Fv'izt)- H',-0 

£> 

(A  2) 

— 1  ridy'  ZzF  Fv  J,z>)  +  E =0  (A3) 

where 

.,  =  i_  p-V 

2((  i/.l„)‘eos<^ 

Ii=p„nA„  cosv? 

C(lCOSl£ 


1)=Kl\Y-\\\ 

/-;=  u(A'l-u 

Solving  equation  (A3)  for  zx  gives 

^-=\F~Gz2-IIFv<lzJ]  i  (A4) 


where 


F~ 

G= 


IF 

if, 

\\\ 

if, 


(1  —  KL)g 


H= 


s 

If, 


Integrating  equation  (A4)  with  respect  to  t  between  tbc 
limits  lT  and  t  and  noting  that  ZiT=Z2T~zr  gives 

zx—zt+Ft— G{zi—  ir)— Fyjzjdr  (Ao) 


where  T=(l—tr). 

Integrating  again  and  noting  that  zlr=z2r=zr  gives 


S|  =  (l  +  (7)(2r  +  Z,t)  + 


Fv'izMr 


dr 


(A6) 
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Substituting  for  i i  anil  _[  in  equation  (A 2)  gives 

1  -C  [(1  +«)(*,+  z,r-22)+  F2~Il££  (^)rfr 

The  motion  of  the  landing  gear  subsequent  to  the  beginning 
of  shock-strut  deflection  is  determined  bv  means  of  a  step- 
by-step  solution  of  equation  (A7).  This  numerical  procedure 
yields  time  histories  of  the  lower-mass  motion  variables 
i2,  and  z2.  from  which  the  motion  variables  for  the  upper 
mass  i:,.  it,  and  zt  can  be  calculated  bv  means  of  equations 
(A4 ) ,  (A5),  and  (A6). 

The  initial  conditions  for  the  step-by-step  procedure  are 


r]|  -) 


1„  =  0—  ) 
'*»()  ~r  / 


where  and  if,  are  the  conditions  of  motion  at  the 

beginning  of  shock-strut  deflection  as  determined  from  the 
solution  for  the  one-degree-of-fr^edom  system. 

Estimated  values  of  the  lo.ve.-mass  velocity  at  the  end  of 
the  Ar  t  time  increment  t  following  the  beginning  of  shock- 
strut  deflection  can  be  obtained  :inm  the  expression 

*?..,=  M—N  (A«) 

or.  as  a  lirst  approximation. 


9*  —7 

~  2  n  =*  1  ~  r 


The  corresponding  displacement  is  given  by 

=  *'.l  =  Zr  +  ±(itl+ir)  (A10) 

After  the  initial  conditions  and  the  conditions  at  the  end 
of  the  first  time  increment  are  established,  a  step-by-step 
calculation  of  the  motion  can  be  obtained  by  routine  opera¬ 
tions  as  indicated  by  the  following  general  procedure  which 
applies  at  any  time  t—ui  after  the  beginning  of  the  process. 
'I'lie  operations  indicated  are  based  on  integration  by  appli¬ 
cation  of  the  trapezoidal  rule: 

— +  (All) 

£  £  2 

-?,7177 +  2  (*2„_|  r«*,)=Ca,_i  +  «i4._,+j  +  5*._2) 


With  the  estimated  values  z* n  and  :*n  the  acceleration  of 
the  lowin'  mass  can  be  determined  bv  substitution  in  the 


appropriate  iutegrodifFcrcntial  equation  for  the  system, 
equation  ( A 7 )  in  the  present  case.  Thus 

=  /(i?  ,  s?  ,  r„)  (A IS) 

In  equation  (A7)  the  integral  expressions  can  also  he  evalu¬ 
ated  by  application  of  the  trapezoidal  rule.  For  example, 
when  FVg  ( z2)  —  m  :/ , 

,AI4’ 

and 

0  Jo  J  0 

•>(Jo  siT,ir+L  zi'dr)  (Ai3’ 

An  improved  value1  for  I  ho  velocity  is  obtained  from  the 
expression 

is  =i»  ~ +  r..  )  (Alt») 

This  value  is  used  in  the  calculation  of  the  estimated  velocity 
ijn+1  and  displacement  z*H . ,  for  the  next  interval. 

If  desired,  improved  values  of  the  displacement  and 
acceleration  for  the  nth  interval  subsequent  to  the  beginning 
of  shock-strut  deflection  can  he  obtained  as  follows: 

=^-,,-1  +  2'  ('%- 1  “ 

+  j  (A  17) 

and 

-**.  =  •>  Pv  .-v  T,)  (A  18) 

where  j<  ,  r2ti,  t„)  is  an  appropriate  equation  for  the  system, 
such  as  equation  (A7). 

With  the  values  of  r2  .  i2  .  and  r2  ,  the  motion  variables  for 
the  upper  mass  .:,n,  i,  .  and  :,u  can  he  calculated  separately 
from  equations  (A4).  ( Ao ) ,  and  (AO),  as  previously  noted. 

In  setting  up  the  numerical  procedure  used  in  obtaining  the 
solutions  presented  in  this  report,  an  evaluation  of  the  errors 
introduced  by  the  procedure  indicated  that  it  would  not  be 
necessary  to  calculate  the  improved  values  of  the  displace¬ 
ment  Zin  (eq.  (Al7))  or  the  acceleration  f2j|  (eq.  (A18)).  How¬ 
ever,  improved  values  of  the  velocity  i2„  were  calculated  by 
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moans  of  equation  (Al6)  for  the  purpose  of  determining 
estimated  values  of  the  velocity  z2  and  the  displacement  z* 
(eqs.  (All)  and  (A12))  for  the  increment,  immediately 
following. 

In  order  to  illustrate  the  application  of  the  method,  a 
tabular  computing  procedure  for  the  solution  of  the  system 
represented  by  equations  (Al),  (A2),  and  (A3)  is  presented  in 
table  I. 

QUADRATIC  PROCEDURE 

In  this  step-by-step  procedure  a  quadratic  variation  of 
displacement  is  assumed  over  successive  equal  finite  time 
intervals  for  the  purpose  of  extrapolating  values  of  the 


motion  variables  from  one  interval  to  the  next.  With  this 
assumption  the  displacement  variation  over  two  successive 
equal  time  intervals  is  completely  determined  by  the  three 
values  of  displacement  at  the  beginning  and  end  of  each  of 
the  two  intervals.  By  writing  the  quadratic  variation  in 
difference  form,  the  velocity  and  acceleration  at  the  midpoint 
of  the  double  interval  can  be  expressed  in  terms  of  the  three 
displacement  values  previously  mentioned.  Substituting  for 
the  velocity  and  acceleration  in  the  differential  equations  for 
the  system  yields  difference  equations  of  motion  in  terms  of 
successive  displacement  values  which  can  be  evaluated 
interval  by  interval. 


TABLE  I 

LINEAR  PROCEDURE 


Row 

Quantity 

Equation 

1 

Procedure f 

© 

© 

1 

| 

+’»-i  +  2  (*!»-j+^s»-i) 

®r 

® 

*=.-,+  2“  (*,-,  +  «.) 

® 

Determined  from  tire  force- 
deflect  ion  charact eristics. 

:© 

Equation  (A  14) 

A  p +  2"  (®  +  ®p] 

® 

JoTfv*(z;")rfr,/r 

E(|iiation  (.415) 

®r  +  -|  [®  +  @p] 

.6) 

-% 

Equation  (.47) 

Given  by  equation  (.47). 

+  2  (*--»-l +  *-'.) 

©p+2'  (®  +  ®rlg 

.  > . 

-J  (*..+«.„) 

+’»  +  {  (ft«-,  +  i»») 

®+ j  I® +©..]« 

(s' 

‘•-n  +  J 

**.  +  2 

©  +  "2  [®+®J* 

!  (1$) 

? 

E<|iiation  (A4) 

Given  by  equation  (.44). 

i  ® 

Zyn 

Equation  (.45) 

(liven  by  equation  (A5). 

'  ® 

Equation  (.46) 

Given  by  equation  (Ati). 

t  Op  denotes 

value  for  previous  time  interval. 
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The  following  derivation  shows  how  the  procedure  can  be 
applied  to  the  determination  of  the  behavior  of  the  landing 
gear  subsequent  to  the  beginning  of  shock-strut  deflection  at 
time  U.  In  order  to  simplify  the  derivation,  internal  friction 
forces  are  again  neglected  in  setting-up  the  equations  of 
motion. 


The  assumption  of  a  quadratic  variation  of  displacement 
with  time  (constant  acceleration)  over  two  successive  inter¬ 
vals,  each  of  duration  «,  permits  expressing  the  velocity  and 
acceleration  at  the  midpoints  of  the  double  interval  (sec 
sketch)  in  terms  of  the  displacement  values  at  the  beginning, 
midpoint,  and  end  of  the  double  interval  by  the  equations 
(see  ref.  5,  p.  16): 


and 


2*+l — 2»-l 
2< 


(AID) 


Sn+l  —  2  Zn  +  Zn-l 


(A20) 


where  i„,  £„,  and  zn  are  the  velocity,  acceleration,  and  displacement  at  the  end  of  the  nth  interval  (r=ne)  after  the 
beginning  of  shock-strut  deformation  and  2„_t  and  2„+1  are  the  displacements  at  the  end  of  intervals  n  —  1  and  n-f  1, 
respectively. 

Substituting  the  difference  relations  for  zh  i2,  zh  and  z2  into  equations  (Al)  and  (A3)  permits  writing  the  equations 
of  motion  for  the  landing  gear  in  difference  form  as  follows: 


and 


77i(2‘n+l— 22l»  +  2l,-l)+|72(2:',+l— ZU-l“22»+l  +  Z*n-\Y  +  B[l~G(ZK~Ztn)]  "+E  =  0 


Sln+l—2Z'n~S'n-l~G(.Z*n+l~2S*n+Z*n-\)~Hf\  Fv,(22n)  +  E] 


where  the  constants  are  as  defined  in  the  previous  section. 
Substituting  for  2in+1  in  equation  (A21)  gives 


■'n+l 


=ir+^riP  [2HVH',-vTHVHygil  Wpn^  +  WSWJ-gA  1H4HV«„+1+5^T,+I)] 


(A21) 

(A22) 


(A23) 


an+i  —  2W  2Z2n—  W  2Z-2n_l  —gt2[F vt(z2n)-\-  A’] 

fl„+i  =  2H  2*?„  +  (W  i  —  ir2)z2n_,  +  2M  vf(^2n)  +  E] 

and 

Equations  (422)  and  (A23)  are  essentially  extrapolation  formulas  which  permit  the  determination  of  values  for  the 
upper-mass  and  lower-mass  displacements  to  come  from  the  values  of  displacement  already  calculated.  These  equations 
thus  permit  step-by-step  calculation  of  the  displacements  as  the  impact  progresses,  starting  with  the  initial  conditions, 
from  which  the  upper-mass  and  lower-mass  velocities  and  accelerations  can  be  determined  by  means  of  equations  (.419) 
and  (A20). 

Since  the  calculation  of  the  displacements  2:  and  ~2  at  any  instant  by  means  of  equations  (A22)  and  (A23)  requires 
values  for  the  displacements  at  two  previous  instants,  the  routine  application  of  these  equations  can  begin  only  at  the 
end  of  the  second  interval  (r=2«)  following  the  beginning  of  shock-strut  deflection.  Before  the  displacements  at  the  end 
of  the  second  interval  can  be  calculated,  however,  it  is  necessary  to  determine  the  displacements  at  the  end  of  the  first 
interval.  These  values  can  be  obtained  from  the  conditions  of  motion  at  the  instant  of  initial  shock-strut  deflection  by 
applying  equations  (A19)  and  (A20)  to  the  instant  t=t,. 
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At  the  instant  of  initial  shock-strut  deflection 


^>».0=S2n.O=^V 


the  foregoing  application  of  the  <lifference  equations  results 
in  identical  values  for  the  upper-mass  displacement  and 
lower-mass  displacement  at  the  end  of  the  first  interval. 
Simultaneous  solution  of  equations  (A25)  gives  the  following 
expression  for  the  displacement  at  the  end  of  the  first  interval: 


-5,-1— —  2r  +  tir  +  -^  5, 


Application  of  the  difference  equations  (A19)  and  (A20) 
to  the  instant  t=tr  (that  is,  n= 0)  gives  the  following 
equations: 


2#  =  I  Z  n  —  l 


\ 

-  _ 1 — S„  =  _i^ 


Since  the  landing  gear  is  considered  as  a  one-degree-of- 
freedom  system  from  initial  contact  up  to  the  instant  (=/T) 


With  the  values  for  z,  and  s„_i,  equations  (A22)  and 
(A23)  permit  the  step-by-step  calculation  of  the  upper-mass 
and  lower-mass  displacements  subsequent  to  the  first  interval 
following  the  beginning  of  shock-strut  deflection.  The 
corresponding  velocities  and  accelerations  of  the  upper  and 
lower  masses  can  be  determined  from  the  calculated  displace¬ 
ments  by  means  of  equations  (A  19)  and  (A20),  as  previously 
noted. 

A  tabular  computing  procedure  illustrating  the  application 
of  the  method  is  presented  in  table  If. 


TABLE  II 


QUADRATIC  PROCEDURE 


Quantity 


Equation 


Procedure  t 


Equation  ( A23) 


Given  by  equation  (A23). 


Equation  (A22) 


Given  !>y  equation  (A22). 


z,  -2  z.  +*.. 


3ft-2  i  +  r 


z'.+.-22'.+2'. 


(l— 2\D+®„ 


t  Op  denotes  value  for  previous  time  interval. 
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Rl'NCE-KUTTA  PROCEDURE 

In  this  step-by-stop  procedure  the  differences  in  the  de¬ 
pendent  variables  over  any  given  interval  of  the  independent 
variable  are  calculated  from  a  definite  set  of  formulas,  the 
same  set  of  formulas  being  used  for  all  increments.  Thus  the 
values  of  the  variables  at  the  end  of  any  given  interval  are 
completely  determined  by  the  values  at  the  end  of  the  pre¬ 
ceding  interval.  Unfortunately,  however,  unless  the  equa¬ 
tions  to  be  integrated  are  relatively  simple,  the  method  can 
become  quite  lengthy. 

The  following  derivation  illustrates  the  application  of  the 
Runge-Kutta  method  to  the  generalized  equations  of  motion 
(eqs.  (21))  for  the  simplified  system  considered  in  the  section 
on  generalized  results.  Since  these  equations  can  be  readily 
reduced  to  the  first  order,  they  can  be  integrated  by  the 
step-by-step  application  of  the  general  equations  given  on 
pages  301  and  302  of  reference  6  for  first-order  simultaneous 
differential  equations. 

The  generalized  equations  for  the  simplified  system  pre¬ 
viously  discussed  (eqs.  (21))  are 

«—«*')•+ tq"=0 

(u/  —  u/)2—u  2=0 

Ui'+u-i— 0 

Inasmuch  as  any  two  of  these  equations  are  sufficient  to 
describe  the  behavior  of  the  system,  only  the  last  two  equa¬ 
tions  are  employed  in  this  procedure.  Those  equations  can 
be  reduced  to  a  first-order  system  by  introducing  the  new 
variable 

w—U\  (A27) 

so  that 

w'=tq"  (A28) 

and  the  equations  of  motion  become 


(w  —  u/f— v2=0 
«•'  -f  w»=  0 


(A  29) 


Solving  equations  (A29)  for  it/  and  «•',  respectively,  gives 
Ui=w—\U3  (A  30) 

(A31) 

Applying  the  general  procedure  presented  in  the  reference 
previously  cited  to  the  simultaneous  equations  (A27),  (A30), 


and  (A31)  gives 

=  —  iq„_,  =  ^-  (k-,  +  2k2+2lc3+k-t) 

Aw=tt\-«^_,  =  J(/, +  24+2/3+/,)  L  (A3  2) 

&Ui  =  ll2n  —  u-“,_i  =  y  (ml  +  2 Hi 2+  2/H3+  Hi,) 

where 

A*i  =  tt’„_iA0 
A‘2=^w„_i+ A0 

A-j =(w„_,+~)  A0 

A-4=(M'b-1  +  4)A0 

AO 

4=-(«2,_,+y)  Af) 

/,=  -(«S„_i+»»j)  A0 


m,=  (tc._1->/w,._l)  A0=tq/_,A0 


Hi,=  [(H'n_l  +4)  —  A0 

With  this  procedure,  iq,  ir,  and  it>  can  be  <  .« -(dated  in 
step-by-step  fashion  from  the  values  for  tin'  preceding  inter¬ 
val,  the  procedure  beginning  with  the  initial  conditions. 
From  these  values,  ?</,  iq",  and  11/  can  be  calculated  by 
means  of  equations  (A27),  (A28),  and  (A30),  respectively. 


APPENDIX  B 


SOURCE  OF  EXPERIMENTAL  DATA 


Following  is  a  brief  description  of  the  apparatus  and  test 
specimen  used  in  obtaining  the  experimental  data  presented 
in  this  report. 

EQUIPMENT 

The  basic  piece  of  equipment  employed  in  the  tests  is  the 
carriage  of  the  Langley  impact  basin  (ref.  7)  which  provides 
means  for  effecting  the  controlled  descent  of  the  test  speci¬ 
men.  In  these  tests  the  impact-basin  carriage  was  used  in 
much  the  same  manner  as  a  conventional  stationary  landing- 
gear  test  jig  (see  ref.  8).  In  order  to  simulate  mechanically 
the  wing  lift  forces  which  sustain  an  airplane  during  landing 
the  pneumatic  cylinder  and  cam  system  incorporated  in  the 
carriage  was  used  to  apply  a  constant  lift  force  to  the  dropping 
mass  and  landing  gear  during  impact.  The  lift  force  in  these 
tests  w  as  equal  to  the  total  dropping  weight  of  2,542  pounds. 

TEST  SPECIMEN 

The  landing  gear  used  in  the  tests  was  originally  designed 
for  a  small  military  training  airplane  having  a  gross  weight  of 
approximately  5,000  pounds.  The  gear  is  of  conventional 
cantilever  construction  and  incorporates  a  standard  type  of 
oleo-pneumatic  shock  strut.  The  wheel  is  fitted  with  a  27- 
inch  type  I  (smooth-contour)  tire,  inflated  to  32  pounds  per 
square  inch.  The  weight  of  the  landing  gear  is  150  pounds. 
The  weight  of  the  lower  mass  (unsprung  weight)  is  131 
pounds. 

In  the  present  investigation  the  gear  was  somewhat  modi¬ 
fied  in  that  the  metering  pin  was  removed  and  the  original 
orifice  plate  was  replaced  with  one  having  a  smaller  orifice 
diameter.  Figure  14  shows  the  internal  arrangement  of  the 
shock  strut  and  presents  details  of  the  orifice.  Other  perti¬ 
nent  dimensions  are  presented  in  table  III.  The  strut  was 
filled  with  specification  AX-VV-0-366B  hydraulic  fluid. 
The  inflation  pressure  with  the  strut  fully  extended  was  43.5 
pounds  per  square  inch.  In  these  tests  the  landing  gear  was 
mounted  with  the  shock-strut  axis  vertical.  Figure  15 
is  a  photograph  of  the  landing  gear  installed  for  testing. 

TABLE  III 

IMPORTANT  CHARACTERISTICS  OF  LANDING  GEAR 
USED  IN  TESTS 

A,,  sq  ft  . .  .  .  0.  05701 

An,  sq  ft - -  -  .  .  .  _  .  _  0.  04708 

A„  sq  ft -  -  -  -  .  0.  0005585 

ro,  cu  ft . .  .  . . .  .  _  0.03545 

JV  lb/sq  ft -  .  _  _  0,264 

h,  ft  . - . .  .  0.5521 

U,  ft  -  - . -  -  - -  2.22604 

H'i,  lb  . .  . . . .  2,411 

W'2,  lb . . .  ...  131 


INSTRUMENTATION 

A  variety  of  time-history  instrumentation  was  used  during 
the  tests.  The  vertical  acceleration  of  the  upper  mass  was 
measured  by  means  of  an  oil-damped  electrical  strain-gag< 
accelerometer  having  a  range  of  ±8gand  a  natural  frequency 
of  85  cycles  per  second.  A  low-frequency  (16.5  cycles  pet 
second)  X  AC  A  air-damped  optical-recording  accelerometer 
having  a  range  of  —  lg  to  6 g,  was  used  as  a  stand-by  instill¬ 
ment  and  as  a  check  against  the  strain-gage  accelerometer 
Another  oil-damped  strain-gage  accelerometer,  having  a 
range  of  ±  12?  and  a  natural  frequency  of  260  cycles  pel 
second,  w-as  used  to  determine  the  vertical  acceleration  of  the 
lower  mass.  The  vertical  displacement  of  the  lower  mass 
(tire  deflection)  and  the  shock-strut  stroke  were  measured 
separately  by  means  of  variable-resistance  slide-wire  poten¬ 
tiometers.  The  vertical  displacement  of  the  upper  mass  was 
determined  by  addition  of  the  strut-stroke  and  tire-deflection 
measurements.  The  vertical  velocity  of  the  landing  gear  at 
the  instant of  ground  contact  was  determined  from  the  output 
of  an  elemental  electromagnetic  voltage  generator.  A  tinu 
history  of  the  vertical  velocity  of  the  upper  mass  was  ob¬ 
tained  by  mechanically  integrating  the  vertical  acceleration 
of  the  upper  mass  subsequent  to  the  instant  of  ground  con 
tact.  Electrical  differentiation  of  the  current  output  of  tin 
strut-stroke  circuit  provided  time-history  measurement* 
of  the  shock-strut  telescoping  velocity.  The  instant  ol 
ground  contact  was  determined  by  means  of  a  micro¬ 
switch,  recessed  into  the  ground  platform,  which  closer 
a  circuit  as  long  as  the  tire  was  in  contact  with  the  platform 

The  electrical  output  of  the  instruments  was  recorded  on  > 
14-channel  oscillograph.  The  galvanometers  were  damper 
to  approximately  0.7  critical  damping  and  had  natural  fre¬ 
quencies  high  enough  to  produce  virtually  uniform  responst 
up  to  frequencies  commensurate  with  those  of  the  measuring 
instrumentation.  A  typical  oscillograph  record  is  shown  it 
figure  16. 

It  is  believed  that  the  measurements  obtained  in  the  test; 


are  accurate  within  the  following  limits: 

Measurement  Accurac; 

Upper-mass  acceleration,  g _ _  0.  ! 

Force  on  upper  mass  lb _ ... -  -  s 50 

Lower-mass  acceleration,  t; -  ±0. : 

Vertical  velocity  at  ground  contact,  fps _ —  ±0. 

Upper-mass  velocity  during  impact,  fps _ _ _  ±0. 

Upper-mass  displacement,  ft _ _ _  ±0.  0 

Lower-mass  displacement,  ft - -  ±0.  0 

Shock-strut  stroke,  ft _  ±0.  0 

Shock-strut  telescoping  velocity,  fps. - - -  ±0. 

Time  after  contact,  sec - ±0-00 
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ANALYSIS  OF  LANDING-GEAR  BEHAVIOR 
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Fioi'rk  14. — Shock  strut  of  landing  gear  tested  at  Langley  impact  basin. 
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F Kil'RK  15. — View  of  landing  gear  and  instrumentation. 
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